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ABSTRACT 

We present results from our spectroscopic study, using the Infrared Spectrograph (IRS) onboard the Spitzer 
Space Telescope, designed to identify massive young stellar objects (YSOs) in the Galactic Center (GC). Our 
sample of 107 YSO candidates was selected based on IRAC colors from the high spatial resolution, high 
sensitivity Spitzer/IRAC images in the Central Molecular Zone (CMZ), which spans the central ~ 300 pc 
region of the Milky Way Galaxy. We obtained IRS spectra over 5 /im to 35 /Ltm using both high- and low- 
resolution IRS modules. We spectroscopically identify massive YSOs by the presence of a 15.4 shoulder 
on the absorption profile of 15 /im CO2 ice, suggestive of CO2 ice mixed with CH3OH ice on grains. This 15.4 
/im shoulder is clearly observed in 16 sources and possibly observed in an additional 19 sources. We show that 
9 massive YSOs also reveal molecular gas-phase absorption from CO2, C2H2, and/or HCN, which traces warm 
and dense gas in YSOs. Our results provide the first spectroscopic census of the massive YSO population in the 
GC. We fit YSO models to the observed spectral energy distributions and find YSO masses of 8-23 M , which 
generally agree with the masses derived from observed radio continuum emission. We find that about 50% of 
photometrically identified YSOs are confirmed with our spectroscopic study. This implies a preliminary star 
formation rate of ~ 0.07 M Q yr~ l at the GC. 

Subject headings: infrared: ISM — ISM: molecules — stars: formation — Galaxy: nucleus 



1. INTRODUCTION 

Our Galactic c enter (GC), at a distance of 7.9 ±0.8 kpc 
dReid et all [2009), is the closest galactic nucleus, observ- 
able at spatial resolutions unapproachable in other galaxies 
(1 pc Rj 26")- The extent of the GC region is defined by 
a region of relatively high density molecular gas («h 2 ~ 10 4 
cm" 3 ; iBally et all 1 19871) . covering the inner 200 pc x 50 pc 
(170' x 40'), called the Central Molecular Zone (CMZ). The 
CMZ produces 5%— 10% of the Galaxy's infrared and Ly- 
man con tinuum luminosity and contains 10% of its molec- 
ular gas dSmith et al.lll978UNishimura et al.lll980tlBallv et all 
1 1 9 871 1 1 98 8h iMorris & Serabynll 19961) . The CMZ exhibits ex- 
treme conditions with high gas temperature, pressure, turbu- 
lence, strong magnetic field strengths and strong tidal shear 
dSerabvn & Morrislll996t iFatuzzo & Melial 120091) . As a re- 
sult, star formation in the CMZ may be altered or suppressed. 

The CMZ, nevertheless, shows several signposts of recent 
massive star formation, such as (compact) H II regions and su- 
pernova remnants. In addition, there are massive young stars 
(ages of -2-7 Myrs; iKrabbe et alJ [19911 iFiger et aljTl999h 
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in three known discrete star clusters - the Central, Quintu- 
plet, and Arches clusters - which make the CMZ distinctly 
different from the Galactic bulge with its p redominantly old 
stellar population dFrogel & Whitfordl ll987l) . Nevertheless, it 
has been unclear how star formation proceeds in this hostile 
environment. There have been several studies in the litera- 
ture that identified young stellar object (YS O) candidates in 
the GC based on in frared photometry (e.g.. iFelli et al.ll2002l 
ISchuller et aill2006t lYusef-Zadeh et al.ll2009l) . The high and 
patchy extinction towards the GC (Ay ~ 30) and its mix of 
young and old stellar populations, however, mean that spec- 
troscopic observations are required to confirm YSO identifica- 
tions. This is because red giants and asymptotic giant branch 
(AGB) stars (also part of the GC stellar population) can look 
like YSOs from bro ad-band photometry, i f they are heavily 
dust attenuated (e.g JSchultheis et alJl2003h . 

The GC provides a unique opportunity to investigate cir- 
cumnuclear star formation with an unprecedented spatial res- 
olution. We announced the first spe ctroscop i c iden tifica- 
tion of massive YSOs in the CMZ dAn et al.l 120091 here- 
after A 09), using the Infrared Spectrograph (IRS .lHouck et alJ 
1 20041) onboard the Spitzer Space Telescope (Werner et al. 
2004). In this paper, we follow up our initial exploration 
of the IRS data set in A09 and refine our methods to iden- 
tify YSOs in the CMZ, aiming at providing a list of spec- 
troscopically confirmed YSOs as tracers of the early stages 
of star formation in the GC. As described and employed in 
A09, our selection criteria for YSOs are based on gas- and 
solid-phase absorption from mid-IR spectroscopy. This in- 
clude s solid-phase a bsorption from the CO2 bending mode 
(e.g.. iGerakines et al]|1999l) a nd gas-phase absorption from 
C7H2, HCN, and CO ? (e.g.. lLahuis & .van Dishoeckl I200TJ: 
iBoonman etail 120031: iKnez et alJ l2009h . We look for sig- 
natures of CO2 ice mixed with a large amount of CH3OH 
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ice. This combination has been o bserved towards high- 
mass YSOs and low-mass YSOs dGerakines et al.1 [l999t 
Ponto ppidan et alj 120081: IZasowski et all 120091: ISeale et all 
20111). but not toward field stars behind molecular clouds 
dGerakines eTaTlll999t iBergin et al.ll200l iKnez etail 120051: 
IWhittet et al.ll2007ll2009h 

In § |2] we summarize the IRS target selection criteria and 
data reduction. In §[3] we describe our spectroscopic identifi- 
cation of YSOs, showing that 15%— 30% of our 107 targets are 
massive YSOs. We measure the extinction for YSOs and pos- 
sible YSOs, along with column densities of solid-phase and 
gas-phase molecular absorbers. In §|4]we examine properties 
of these YSOs and possible YSOs, and derive a preliminary 
estimate of the star formation rate in the GC. 

2. METHODS 

In this section we describe procedures for the sample selec- 
tion, spectroscopic follow-up observations, and IRS data re- 
duction. Parallel information on these subjects can be found 
in A09, but here we repeat this for the reader's convenience 
with additional details where there has been improvements in 
the data reduction steps. 

2.1. Spitzer/IRS Sample 

Our 107 spectroscopic targets (Table [2} we re sele cted from 
the GC point source catalog (Ramirez et al. 2008) extracted 
from the IRAC images of the CMZ dStolovv et al.l2006]) made 
using the Infrared Array Camera (IRAC: iFazio et aiTl2004b 
onboard the Spitzer Space Telescope. These images cover 
280 pc x 200 pc in the four IRAC channels (3.6 pm, 4.5 pm, 
5.8 pm, and 8.0 pm) with uniform high sensitivity. Compared 
to earlier imaging surveys of this regio n, such as that fr om 
the Midco urse Space Experim ent (MSX; iPrice et aTll2001l) or 
ISOGAL dOmont et al.l2003l) . the IRAC images have a higher 
spatial resolution (« 2" vs. > 6" of earlier surveys), which 
has led not only to a better estimate of source fluxes, but also 
to more accurate source positions for follow-up spectroscopic 
observations. 

The spectroscopic samp le was selected using IRAC color 
criteria based on the Whi tney et al.l d2004l) study of the giant 
H II region RCW 49. 1 Whitney et al.l determined the locations 
of YSOs with 2.5 M Q , 3.8 M , and 5.9 M Q on the IRAC 
color-mag nitude diagrams, usin g radiative transfer models de- 
scribed in Whitnev et al. (2003J). From this we chose an ini - 
tial color criterion ([3.6]- [8.0] > 2.0: IWhitnev et al.ll2004l) . 
We added a latitudinal constraint (\b\ < 15') to increase the 
probability that the objects a re located at the distance of the 
GC (8 kpc; iReid et al J 120091) rather than in one of the sev- 
eral intervening spiral arms along the line of sight. We note 
that the range of this latitude selection is about 5 times larger 
th an the scale height of ph otometric YSO candidates (~ 7 pc) 
in lYusef-Zadeh et al.1 d2009l) . These color and position con- 
straints provided an initial sample of 1207 objects from the 
GC point source catalog. 

We combined the IRAC photometry wit h Two Micron 
All S ky Survey (2MASS) photometry (7fflf,; ISkrutskie et al l 
2006) and ISOGAL 7 pm and 15 pm point-source cata- 
logs. Note that 24 /im Multiband Imaging Photometer for 
Spitzer (MIPS) observations at the GC were not available at 
the time when our IRS sample was chosen. Among the ini- 
tially selected 1207 objects, 336 had photometry in at least 
5 bandpasses, which allowed reliable spectral energy distri- 
bution ( SEP) fitting. We used the SED fitting tool devel- 
oped by iRobitaille et al.l d2007h . which makes use of a grid 



of 200,000 YSO models dRobitaille et all 120061) to estimate 
Y SO parameters, such as the mass of a central object (see 
§ 14. lb . For those objects with photometry in fewer than 
five bandpasses , we instead applied color constraints based 
on the work of IWhitney etail d200l of [3.6] -[4.5] > 0.5, 
[4.5] - [5.8] > 0.5, and [5.8] - [8.0] > 1 .0. 

The SED fitting and additional IRAC color constraints nar- 
rowed our sample down to about 200 objects, which were then 
further inspected using IRAC three-color images. The sources 
were evaluated by their distinctiveness (i.e., whether sources 
are easily distinguishable from the background) and their lo- 
cal background emission. Among 200 objects examined, 1 12 
were found to exhibit the necessary distinctiveness within the 
IRS slit widths. 

A literature search was performed on the 112 objects, yield- 
ing matches to 43 previously studied sources: 25 sources were 
previously-identified photometric YSO candidates, 4 were 
OH/IR stars, one was a Wolf-Rayet star, and the remaining 

13 sources were others (e.g., radio sources, X-ray sources, 
etc). Note that nearly 60% (« 25/43) of the objects had been 
selected as YSO candidates by other methods. The 4 OH/IR 
stars and the Wolf-Rayet star were discarded from the final 
sample, giving a total of 107 massive YSO candidates. The 
spatial location of the 107 massive YSO candidates of our 
sample is shown in Figure Q] 

In spite of our efforts to exclude OH/IR stars from our YSO 
sample, we later realized that several of our targets appear to 
be coincident with a stellar maser source and/or a long-period 
variable (SSTGC 284291, 425399, 564417, 619964, 660708, 
696367, and 711462 ). SSTGC 517724 is n ow identified as 
an OB supergiant by Mauerh an et al.1 (12010b . These sources 
known not to be YSOs have been helpful in refining our spec- 
troscopic YSO selection criteria. 

Massive YSOs are our pri mary targets for the follow-up 
spectroscopic observations (§ 12.2) because our adopted color 
selection criteria set the lower lim it on the mass of the cen - 
tral object to be M* > 2.5 M Q dWhitnev et alj|2003l 120041) . 
In addition, the source confusion li mit in the input GC point 
source catalog (see Figure 12 in Ra mirez et al.l 120081). to- 
gether with the 8 kpc distance to the GC dReid et al.ll20"09h 
and Ay ~ 30 mag of visual extinction, limits us to detecting 
YSOs with masses > 6 M (see §0). Thus any YSO we iden- 
tify in this paper is a massive YSO. 

2.2. Observations 

Our IRS observations with a total integration time of 
56 hours were carried out in May and October 2008 
(see Table |2} as part of Spitzer Cycle 4 (Program ID: 
40230, PI: S. Ramirez). We observed our 107 targets 
with both high- and low-resolution IRS modules: short- 
high (SH; 9.9 /mi-19.6 pm, A/AA - 600), long-high 
(LH; 18.7 /im-37.2 pm, A/AA ~ 600), short-low (SL; 
5.2 pm-U.5 pm, A/AA ~ 60- 127), and long-low (LL; 

14 pm-38 pm, A/AA - 57- 126). 

In Table [3] we list IRS modules used in the current analysis 
for each of the sources that are spectrosco pical ly identified as 
a YSO or possible YSO in this paper (see § 13. It . Most of these 
targets were observed with all of the four IRS modules. We 
did not obtain spectra with the SL module for some of the 107 
YSO candidates, including possible YSO SSTGC 610642, be- 
cause of saturation in the IRS peak-up arrays (see below). For 
a few sources we rejected data in the first order of LL (LL1; 
19.5 /im - 38.0 pm) because a large fraction of pixel values 
were flagged as invalid. 
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FIG. 1. — Spatial distribution of 107 IRS targets on the IRAC 8.0 /im image i Stolovy et al. 2006). The image shows the entire CMZ covering approximately 
100' X 40' centered on the GC. Our IRS targets (shown in circles) were selected from the point sources of this survey and they are uniformly distributed over the 
CMZ. 



We divided our sample into four subsamples according to 
their IRAC [8.0] magnitudes: [8.0] < 6 mag (N = 30 objects), 
6 mag < [8.0] < 7 mag (N = 28), 7 mag < [8.0] < 8 mag 
(N = 28), and [8.0] > 8 mag (N = 21). Exposure times were 
determined for each brightness subsample to achieve a signal- 
to-noise (S/N) ratio of at least 50 in SH and SL, and a mini- 
mum S/N of 10 in LH and LL. Our exposure times are 6 sec- 
120 sec in SH, 6 sec-60 sec in LH, 6 sec-14 sec in SL, and 
6 sec in LL modules. Subsamples were further grouped based 
on spatial location. These groupings allowed us to observe 
107 sources using nine "fixed cluster" target observations, a 
strategy which proved to greatly increase the observing ef- 
ficiency by reducing overheads due to telescope movement. 
Each object was observed in the IRS staring mode with 4 ex- 
posures per source (2 cycles) to properly correct for bad pix- 
els. 

Our observations were carried out without specific IRS 
peak-up sequen ces, since target c oordinates were accurate 
enough (< 1"; Ram irez et al.l l2008) for our science goals. In 
addition, the background at the GC is too high for a peak-up 
sequence to work even using a 2MASS source. In fact, the 
background at the GC is so high that constraints were placed 
on the observing dates (which determine the telescope roll an- 
gle) to avoid saturation of the IRS peak-up arrays. Such satu- 
ration leads to incorrect droop corrections in the standard IRS 
pipelines and causes various defects on SL frames (where the 
peak-up arrays are located). The SL observations were not 
carried out for 21 out of 107 targets because of saturation in 
the peak-up arrays regardless of the date of observation. 

Multiple off-source measurements at several different loca- 
tions were carried out to derive background spectra around 
each target, because strong and spatially variable background 
at the GC can affect resulting line and/or continuum emis- 
sion from the source. Since the high-resolution slits are not 
long enough to take both source and background measure- 
ments simultaneously, we located four background positions 
around each target (~ ±1 ' offsets in right ascension, ~ ±1' 
offsets in declination). Specific background positions for both 
SH and LH were determined to avoid background sources and 
to properly interpolate background emission near the source 
position over a ~ 1' scale. The longer slit sizes of the low- 



resolution modules permit background measurements along 
the on-source slit; we also identified two additional back- 
ground positions that are ~ ± 1 ' away in the direction perpen- 
dicular to the SL or LL slit. These dedicated background slits 
for SL and LL were centered on two of the high-resolution 
background positions. 

Figure |2] displays the IRS slit positions for on-source (left) 
and off-source (right) measurements for one of our sources 
(SSTGC 797384). Source and background spectra were taken 
consecutively to minimize zodiacal light and instrumental 
variations. Each order of SL or LL was used to observe a tar- 
get, and different orders cover different parts of the sky near 
each target. The low ecliptic latitude of the GC restricts the 
LL slits to a position angle (PA) of ~ ±90° and the SL slits 
to a PA of- 0° or 180°. 

2.3. Data Reduction 

We began reducing the high-resolution IRS spectra from 
the basic calibrated data (BCD), while we started with coad- 
ded products (post-BCD) for the low-resolution spectra. We 
used the S18.7 version of the IRS pipeline for both. On the LH 
frames we applied the DARKSETTLE0 software package to 
correct for non-uniform dark currents. We corrected for rog ue 
pixel values using the SSC software package IRSCLEAN 1 . 
We only applied campaign rogue masks (Campaigns 50 and 
55 for the spring and autumn runs, respectively), except in SL, 
where we applied our own edited version to mask out hot pix- 
els at — 10 fim in addition to campaign rogue pixels. We then 
used SPICE 1 to extract target and background spectra, and 
further corrected high-resolution spectra (SH, LH) for fringe 
patterns using the IRSFRINGE 1 package. 

Four background spectra per target were extracted for the 
high-resolution (SH and LH) observations as four off-source 
pointings were obtained per target. For the low-resolution (SL 
and LL) observations, two background spectra were extracted 
from observations at the same positions as the high-resolution 
background observations and two background spectra were 
extracted along the on-source slit observations. For some ob- 

10 The SSC software packages can be found at 

http://ssc.spitzer.caltech. edu/dataanalysistools/tools/ 
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FIG. 2. — IRAC 8.0 (im image showing a 4' X 4' field of view, centered on one of our IRS targets (SSTGC 797384). Left: positions of all IRS slits (SL: red, LL: 
yellow, SH: magenta, and LH: green) for on-source measurements. Right: Off-source measurements, showing background positions around the source. Similar 
slit formations were adopted for all of the spectroscopic targets. The four off-source pointings were observed to derive a background spectrum for each source, 
because of strong and spatially variable background towards the GC. 
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FIG. 3. — IRS spectra of spectroscopically identified YSOs. Orange lines are low-resolution (SL, LL modules) spectra, and green lines are high-resolution (SH, 
LH modules) spectra. The high-resolution spectra were scaled to match the flux in low-resolution modules (see text). 



jects, the high-resolution background slits were not coinci- 
dent with the on-source low-resolution slits. For these, we 
inspected the slit positions on an IRAC/MIPS composite im- 
age to determine the extraction position along the on-source 
SL or LL slit that was closest in flux and position to the high- 
resolution background slits. 

For each module and for each target, we were able to ex- 
tract four background spectra, that we used to estimate the tar- 
get background by making a linear interpolation of the back- 



ground flux at the source position at each wavelength. Our 
interpolation scheme estimates the background flux gradient 
over a ~ 1' angular scale, since each background pointing is 
~ 1' away from the science target. If the background emis- 
sion is varying over a smaller angular scale, then background 
subtraction would be more uncertain. 

We found that spectra extracted from various IRS modules 
usually do not match with each other at overlapping wave- 
lengths, primarily due to the different sizes of the slit en- 
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FIG. 4. — Cont'd. Same as in Figuref5] 



trances. The SL and LL modules have 3.7" and 10.7" slit 
widths, respectively, while the slit entrances of SH and LH 
modules are 4.7" x 11.3" and 11.1" x 22.3", respectively. 
Therefore, contamination by point sources and/or extended 
emission at the GC can easily lead to a flux mismatch among 
the various IRS modules. 

To obtain internally consistent fluxes from all of the IRS 
modules, we scaled the spectra to match the fluxes from the 
second order of LL (LL2; 14 /im - 21 /im). On the longer 
wavelength side, we scaled the LL1 spectrum to the LL2 
spectrum by estimating a median flux ratio for the two mod- 
ules in the overlapping wavelength region. We masked known 
emission features and rejected points that were more than 3 a 
away from the median flux ratio. On the shorter wavelength 
side, the scaling was done in a step-by-step fashion. We first 
scaled the SH spectrum to the base flux of the LL2 spectrum 
(before correcting for order tilts; see below). We then matched 
the spectrum in the first order of SL (SL1; 7.4 /im - 14.5 fim) 
to the SH spectrum, then scaled the flux in SL3 (bonus order; 
7.3 /xm-8.7 /jm) to the SL1 spectrum, and finally scaled the 
flux in the second order of SL (SL2; 5.2 //m - 7.7 /im) to SL3. 
We can describe this concisely as LL2 — > LL1, and LL2 — > 
SH — > SL1 — > SL3 — > SL2. If a source was not clearly sep- 
arated from extended background emission in LL2, we opted 
to choose the SL1 spectrum as the base flux for scaling. In 
this case, the flux calibration was done in the following se- 
quences: SL1 -> SH ->■ LL2 -> LL1, and SL1 -> SL3 -> 
SL2. 

Scaling factors applied to each module are shown in Ta- 
ble [3] The second column in Table [3] shows which mod- 



ule was used as the baseline for the flux calibration for each 
target. Either large or small scaling factors are found in 
LL for SSTGC 360559, SL for SSTGC 773985, and LL for 
SSTGC 801865. These objects are faint ([8.0]= 7.6 mag and 
8.8 mag for SSTGC 360559 and 801865, respectively, while 
SSTGC 773985 was not detected in this bandpass) on top of 
bright or saturated background emission on MIPS [24] im- 
ages. As a result, their mid-IR fluxes are heavily contami- 
nated by background emissions at A > 20 /im, leading to an 
over-estimation of flux from the target slit. 

Order-tilt features remained in about 30 high-resolution 
spectra, after applying IRSFRINGE to SH and LH and 
DARKSETTLE to LH. To remove this artifact, we applied a 
1 st order polynomial to each high-resolution spectral order to 
force it to match the re-scaled low-resolution spectra. In ad- 
dition, three sources showed scalloping features in their high- 
resolution spectra; a 2 nd order polynomial was applied to cor- 
rect for this artifact. Individual spectra from various orders 
were then merged together using a linear ramp. 

Some data were excluded from the analysis due to problems 
in a particular spectrum, such as saturation, excess bad pixels, 
or a poor match in background level. The last column in Ta- 
ble[3]lists any excluded data. Each target was observed at two 
different nod positions in each module; any nod positions that 
were excluded for a particular module are also given in the last 
column of Table [3] A cardinal point (NSEW) given in paren- 
theses in this column indicates that the background spectrum 
offset in that direction from the source was not included in the 
background determination for that module. 

Figures [3}{7]display spectra resulting from the above proce- 
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FIG. 5. — IRS spectra of possible YSOs. Line colors are the same as in FigurefS] 



dures; orange lines are low-resolution spectra, and green lines 
represent high-resolution spectra. Only sources we spectro- 
scopically identify as a YSO or possible YSO are shown in 
Figures [3}|6] Spectra of known OH/IR, long-period variable, 
or OB supergiant stars in our sample are shown in Figure [7] 
for comparison. 

There are several sources of flux uncertainty in our spec- 
tra: statistical, calibration (difference between different nods), 
and the varying background (this last is usually largest). This 
can cause spectral features observed in emission in the back- 
ground appear in absorption in some spectra, such as the 
1 1.3 /urn polycyclic aromatic hydrocarbon (PAH) feature (e.g. 
SSTGC 304239), H 2 emission at 17.0 pm (e.g. SSTGC 
761771), or forbidden lines such as 12.8 /im [Ne II] or 18.7 
/zm [S III] (e.g. SSTGC 670953). We have estimated the un- 
certainty due to background subtraction by comparing results 
derived by excluding one of the four background pointings 
from the interpolated background spectrum. Throughout our 
analysis, we have added these uncertainties in quadrature to 
derive final uncertainties in measured quantities such as the 
CO2 ice column density. 

3. RESULTS 

3.1. Spectroscopic Identification of Massive YSOs in the GC 



To study the 15 /im CO2 ice absorption profile, we fitted 
five laboratory spectral components to the feature in all our 
107 targets, following the same procedure described in A09. 
FigureslMTOlshow the CO2 ice decomposition for our YSOs, 
possible YSOs, and known stars for comparison, respectively. 
We describe below how we selected YSOs and possible YSOs 
based on this procedure. 

First, we set a local continuum over 14.3 pm < A < 
16.5 pm using a 3 rd order polynomial to derive the optical 
depth. Then we used the mode ling technique and laboratory 
data in Pontoppidan et al. (2008) to decompose the absorption 
profile with five laboratory spectral components; these are po- 
lar C0 2 (C0 2 :H 2 = 14 : 100 at 10 K; dotted line, centered at 
- 15.3 Aim), apolar C0 2 (CO:C0 2 = 100 : 70 at 10 K; dotted 
line, centered at ~ 15.1 fim), pure CO2 (15 K; blue shaded), 
diluted C0 2 (CO:C0 2 = 100 : 4 at 10 K; black solid line), 
and 15.4 A*m shoulder CO2 (modeled with two Gaussians in 
wavenumber space; orange shaded). We found a best-fitting 
set of m odels from the no n-linear least squares fitting routine 
MPFIT dMarkwardtll2009l) . 

Fitting results are shown in Table |4] The CO2 ice 
column densities were estimated from the integrated ab- 
sorption, adopting the integrated line strength A = 1.1 x 
10" 17 cm molecule -1 dGerakines et al.l 19951) . Background un- 
certainties were estimated by creating spectra with one of four 
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FIG. 6. — Cont'd. Same as in Figure|5] 



background positions excluded from the interpolated back- 
ground spectrum and then comparing the column den- 
sities derived from these spectra. We added these uncertain- 
ties in quadrature to the uncertainties from comparing column 
densities from spectra at the two nod positions, and to uncer- 
tainties in column densities due to the statistical uncertainties. 
The Xtot an d Mot in Table |4] represent the total chi-square of 
the fit and the number of data points used in this fit. The 
goodness of fit is generally poor, implying either underesti- 
mated flux errors or our lack of knowledge of individual CO2 
ice models. Nevertheless, the ice decomposition still provides 
useful information on the nature of YSOs, as shown below. 
For comparison, fitting results for some known stellar sources 
are included in Table |4] 

Our primary method of identifying YSOs from our IRS ob- 
servations is the CO2 ice absorption profile at 15 /zm, which 
is obse rved to have a different spectral shape in and around 
YSOs (lEhrenfreund et al.lfl999l IDartois et al.lll999al) . High- 
spectral reso lution observations of many massive YSOs in 
our Galaxy (Gerakines et al. 1999) and in the Large Mag- 
ellanic Cloud dSeale et al.ll201 ll) found a "shoulder" at 15.4 
/zm on the CO2 ice absorption profile. This 15.4 /zm 
shoulder is thought t o be due to the presenc e of CH3OH- 
rich C O? ice grains (Ehrenfr eund et all 119991: IDartois et al. 
Il999al) . Detailed fitting of the 15 /an CO2 ice profile shows 



that the 15.4 /zm shoulder is weaker in low- mass protostars 
(Pontoppidan et al. 2008; Zasowski et al. 2009) and is not de- 
tected towards field s t ars behind several molecular clouds 
dGerakines et al.lll999t iBergin et all 120051: iKnez elaTI l2005h 
IWhittet et alJ 120071 120091) . Analysis of the 15 /zm C0 2 ice 
profile along the lines of sight to the Central Cluster and to 
two dusty WC9 stars in the Quintuplet Cluster demonstrates 
that none of these three GC spectra shows a 15.4 /zm shoul- 
der on the 15 /im CO2 ice absorption profile (Gerakines et al. 
[19991) . Thus, the presence or absence of the 15.4 /im shoulder 
is an empirical - and quantitative - way in the GC of distin- 
guishing YSOs from AGB stars behind molecular clouds. 

The 15 /im CO2 ice absorption profiles displayed in Fig- 
ures [8}j9] for 35 of our YSOs or possible YSOs (see below) 
show two absorption peaks, at 15.15 /im and 15.4 /im. Many 
previously studied YSOs show a double-peaked absorption 
profile, but with pe aks at shorter wavelengths of 15.10 /im 
and 1 5.25 /zm ("e.g..lGerakines et al.lll999t lp"ontorjpida n et al.1 
l2008tlSeale etalJl2011l) . Double-peaked absorption at 15.10 
/im and 15.25 /zm is ascribed to pure CO2 ices result- 
ing from the crystallization of heated H2Q-r ich ices (e.g., 
Gerak ines et al.lll999t iPontoppidan et al.ll2008l). By contrast, 
Ehrenfreund et al.1 (1 1999b and IDartois et all (11999a!) interpret 
CO2 ice absorption peaking at 15.15 /zm as due to CO-rich 
CO2 ices and absorption peaking at 15.4 /zm as arising in 
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FIG. 7. — IRS spectra of known stars (non-YSOs) in our target sample. SSTGC 517724 is an OB supergiant star (Mauerhan et al. 2010) and the other targets 
are OH/IR stars or long period variables. Line colors are the same as in Figure[5] 



CH 3 OH-rich C0 2 ices. 

We selected YSOs by requiring that the model fit to the ob- 
served 15 /im CO2 ice profile significantly improves when the 
15.4 /im shoulder is included in the model. We calculated the 
reduced \~ for fitting a four-component model (excluding the 
15.4 /im shoulder) to the 15 urn CO2 profile, as shown in the 
bottom panels in Figures |8]-|9| We then calculated the reduced 
X 2 for fitting the five-component model (including the 15.4 
/mi shoulder). Finally, we calculated Ax 2 , equal to the re- 
duced x 2 for the four component model minus the reduced 
X 2 for the five component model. We also required that the 
optical depth from the 15.4 /im feature be more than 0.05, a 
limit set by the IRS flat field uncertainty. This corresponds 
to a column density for the 15.4 /im CO2 ice component of 
^(shoulder) » 0.5 x 10 17 cm" 2 . 

We illustrate our YSO selection in Figure [TT] where we plot 
Ax 2 vs. Af co i(shoulder). We conclude that a GC source is a 
YSO if A X 2 > 2 and AU (shoulder) > 0.5 x 10 17 cm" 2 . We 
define a GC source as a possible YSO if < Ax 2 < 2 and 
/Vco^shoulder) > 0.5 x 10 17 cm" 2 . We visually inspected pos- 
sible YSO spectra, and excluded some spectra as clearly non- 
YSO: these are SSTGC 440424 (weak 15 /im C0 2 absorp- 
tion), SSTGC 564417 (OH/IR star), SSTGC 619964 (vari- 
able star), SSTGC 696367 (OH/IR star), SSTGC 660708 
(OH/IR star), SSTGC 732531 (15 /im C0 2 absorption not sig- 
nificantly different between source and background spectra), 
SSTGC 738126 (weak 15 /im C0 2 absorption). We consider 
all other GC sources not to be YSOs. These cutoff values of 
Ax 2 and A^ co i(shoulder) closely agree with the YSO classifi- 



cation that three of us (DA, SR, KS) did by visually inspecting 
the IRS spectra of all 107 targets. 

Our spectroscopic classification of the 107 GC targets is 
shown in the fifth column of Table [2] We conclude that 16 
sources are YSOs ("yes" in me fifth column of Table |2]i and 
19 sources are possible YSOs ("maybe" in the fifth column 
of Table 01. The remaining columns in Table |2] show cross- 
identifications of our IRS sample w it h earl ie r photometry- 
based YSO selections in iFelli et a l. (2002), Schull er et all 
(2006), and lYusef-Zadeh et al.l (120091) . We describe these 
cross-identifications in § 14.41 

The strength of the 15.4 /im peak in our sources is simi- 
lar to that of the well- studied embedded massive YSO W33A 
dGerakines et al.|[T999l) . It is ascribed to a Lewis acid-base in- 
teracti on of C0 2 (the Lewis acid) with CH3OH (Da rtois et al.l 
1999a). Other species could be acting as a base as well, but 
CH3OH is preferred due to its high abun dance toward W33A , 
which is 5%-22% relative to solid H 2 dDartois et al.ll999bl) . 
Two other massive YSOs (AFGL 7009S, AFGL 2136) show a 
prominent 15.4 /im peak, an d indeed these sources have high 
CH3O H abundances as well dDartois et al.ll 1 999bt iGibb et all 
120041) . We therefore suggest that the GC YSOs and pos- 
sible YSOs may also have high solid CH3OH abundances. 
Although the origin of the large quantities of CH3OH in 
the pr eviously studied mas sive YSOs is not fully under- 
stood (Dartois et al.l Il999al) . all lines of sight with high 
solid CH3OH abundances are associated with star formation, 
strengthening the argument that the sources studied in this pa- 
per are indeed YSOs. 
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3.2. Gas-phase Absorption 

Many of our YSOs have gas-phase absorption from C2H2 
(13.71 fjxn, v 5 = 1-0), HCN (14.05 fim, u 2 = 1-0), and/or 
CO2 (14.97 fim, 1/2=1 -0). These gaseous bandheads have 
been detected t oward other massive YSOs, tracing warm and 
dense gas (e.g.jLahuis & van Dishoeckl2000T:lBoonman et alj 
120031: fenez et alj|2009l) All GC sources with these gas ab- 
sorption bands have been already identified as YSOs through 
the strength of the 15.4 fim CO2 ice shoulder, thus strength- 
ening our identification. An AGB star can show either CO2 
gas or C2H2 gas but not both, because while O-rich AGB 
stars sometimes show CO2 gas in emission or absorption 
dJusttanont et al. 1998), C 2H2 gas a bsorption is found only in 
C-rich AGB stars (lAoki etafll 19991) . 

Figure Q~2] shows relative intensity spectra for the nine 
YSOs in our sample that show gas-phase absorption from 
at least one of these species. Three YSOs (SSTGC 524665, 
SSTGC 797384, SSTGC 803 187) presented in A09 are shown 
in Figure [T2l together with six additional YSOs with signifi- 
cant gas-phase absorption. The relative intensity was deter- 
mined by using a second order polynomial to set a local con- 
tinuum at 13.30 ^m < A < 14.55 ^m for C 2 H 2 and HCN, and 



at 14.77 fim < A < 15.06 fim for C0 2 . 

As in A09, we used model spectra from Cami et al. 
(1201 01). These model s are based on the HITRAN04 linefist 
dRothman etail 120051) for C 2 H 2 and HCN, and based on 
HITEMP dRothman etaill 19961) for C0 2 . We did not include 
isotopes in the computation because of the limited parame- 
ter span in the model grids. However, even a relatively high 
isotopic fraction in the GC ( 12 C/ I3 C - 25; iWannieri [19801 
iGtisten et al.lll985l) has a negligible impact on the model fit- 
ting. Best-fitting model values of the excitation temperature, 
T ex , and the gas-phase column density, N CO i, were found by 
searching for the minimum \ 2 of the fits over 100 K < T ex < 
1000 K in steps of Ar ex = 100 K, and 15 < logJV co i < 18 for 
C 2 H 2 , 16 < logJV col < 18 for HCN, and 16 < logJV co i < 22 
for CO2 with intervals of 0.1 dex. Errors in these parameters 
were estimated from A^ 2 , where 1 a measurement errors were 
taken from the scatter in flux. Systematic errors from back- 
ground subtraction and nodding differences were then added 
in quadrature. We tested varying covering factors (the fraction 
of the background continuum source covered by the compo- 
nent in question), but found that the best-fitting value is equal 
to or close to unity. We adopted a Doppler parameter of 3 km 
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Figure [12] shows the best-fitting models for each molecular 
species in red lines. Individual absorption lines are marked 
with vertical bars if they were identified by three of us (DA, 
SR, KS) by visually inspecting the IRS spectra of all 107 
targets, independent of the model fitting. Some of the lines 
were marked undetected (e.g., C 2 H 2 of SSTGC 761771) be- 
cause of a low signal-to-noise ratio of its spectrum. The best- 
fitting model excitation temperatures (r ex ) and column densi- 
ties (A'coi) of identified lines are listed in Table [5] All objects 
with identified gas-phase absorptions are selected as YSOs 
through the detection of the 15.4 fim shoulder component of 
the 15 /im CO2 ice absorption feature. 

3.3. Extinction 

The extinction for our sources can be derived from the op- 
tical depths of 9.7 /im and 18 /im silicate absorption features 
in the IRS spectra. We derived two estimates of the dust ex- 
tinction: one using the low -resolution modules SL+LL [here- 
after Ay (SL+LL)], and one using the high-resolution mod- 
ules SH+LH [hereafter A v (SH+LH)]. The determination of 
Ay (SL+LL) takes both the 9.7 /im and 18 /im silicate fea- 
tures into account. The high-resolution data do not include the 
short wavelength side of the 9.7 /im silicate feature, and so Ay 
(SH+LH) is mainly constrained by the 18 /im silicate feature. 
The 18 /im feature is broader and shallower than the 9.7 /im 
absorption, so it provides a weaker constraint on Ay. Ay 
(SH+LH), however, provides a useful diagnostic when SL is 
not available, as many sources near Sgr A do not have SL data 
due to saturation in the peak-up arrays (e.g., SSTGC 610642). 
Since the high-resolution spectra were scaled to the flux in the 



low -resolution modules, Ay (SL+LL) and Ay (SH+LH) are not 
independent from each other. 

To determine the dust extinction, we model the 5 /im - 
32 fim spectrum by simultaneously fitting the underlying con- 
tinuum, the silicate dust features centered at 9.7 /im and 18 
/im, and the 13 /im librational H2O ice absorption (see Fig- 
ure 4 in A09). The entire silicate extinction curve, derived us- 
ing the GCS 3 spectrum from the Infrared Spac e Observatory 
(ISO) Short Wavelength Spectrometer (SWS) (iKemper et al.l 
2004), is characterized by the optical depth at 9.7 fim, T9.7. 
We adopted the la boratory spectrum of pure amorphous H2O 
ice at T = 10 K (Hudgi ns et alJll993h to model the 13 fim 
librational H2O absorption. This shallow absorption is not 
well-constrained, however, so that the resulting column den- 
sity of H2O ice, iV co i(13 fim), is uncertain. We used a second- 
order polynomial to simulate the overall shape of the SED 
plus grey extinction in the absence of silicate and H2O ab- 
sorption. Before p erforming this non-linear least squares fit 
(Markwardt 2009) we masked molecular absorption features 
at 5.5 fim < A < 7.5 fim, PAH emission at ~ 11.3 /im, CO2 
ice absorption at ~ 15 fim, strong emission lines, as well as 
the noisy bottom part of the 9.7 fim silicate feature (9.3 fim 
< A < 10.1 fim). We derived Ay (SL+L L) from T9.7 by adopt- 
ing Ay/797 = 9 (Roche & Aitken 1985), the value measured 
for lines of sight towards the GC. We derived Ay (SH+LH) in 
the same way, except that we modeled the 10 fim - 32 fim 
high-resolution spectra instead. 

The uncertainties in Ay (SL+LL) and Ay (SH+LH) are dom- 
inated by the uncertainty in choosing the continuum. We es- 
timated these uncertainties by comparing results where the 
continuum was derived from the same wavelength regions in 
all spectra to results where each continuum was set interac- 
tively. Applying a second-order polynomial for a continuum 
generally results in a good fit over 10 fim < A < 32 fim, but 
underestimates fluxes at < 8 fim, which may be due to under- 
subtraction of background PAH emission at ~ 7 .7 fim. We 
followed the prescription in Booge rtet al.l (12008) to force the 
continuum to match (by eye) the observed flux at ^5.5 fim 
and ~ 7.5 fim and to set an approximate flux at ~ 30 fim. 
We tried a number of interactive continuum settings, but this 
approach generally results in a much worse agreement of the 
model fits with observed flux over 10 fim < A < 32 fim. We 
took this as an upper ler boundary of Ay (SL+LL). Errors in 
Ay (SH+LH) include statistical uncertainties, where we took 
the scatter of points at 20 fim < A < 30 fim with respect to 
a second-order polynomial as the flux errors over the entire 
wavelength range, added in quadrature to uncertainties from 
varying the background subtraction and uncertainties between 
the two nod positions. 

Table|6]shows Ay (SL+LL) and Ay (SH+LH) estimates for 
YSOs and possible YSOs. We compare these two extinction 
estimates in the upper panel of Figure Q~3] The extinction for 
each YSO is a combination of extinction along the line of 
sight to the GC and extinction intrinsic to the YSO; sources 
which are not YSOs will not always have intrinsic extinc- 
tion. There is a good correlation between Ay (SL+LL) and 
Ay (SH+LH) as illustrated in Figure [13] This is expected, 
because the high-resolution spectra are scaled to the low- 
resolution spectra, and so the two methods are not completely 
independent from each other. The weighted mean difference 
is (Ay(SH + LH)- Ay (SL + LL)) =+0.38 ±0.65 mag for both 
YSOs and possible YSOs. The (unweighted) rms difference 
is 1 1 mag, compared to the formal uncertainties of ~ 9 mag 
from both axes. 
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are selected as YSOs through the detection of the 15.4 fim shoulder component of the 15 fim C0 2 ice absorption feature. 



The lower panel in Figure Qj] shows a comparison of Ay 
(S L+LL) to Ay (foregroun d) derived from the extinction map 
in lSchultheis et al.l d2009t) . The latter is based on the 2MASS 
and IRAC color-magnitude diagrams of GC red giant branch 
stars within 2' from each source. The errors are the rms differ- 
ence of Ay (foreground) derived at the positions of four back- 
ground pointings. As seen in the figure, Ay (SL+LL) is sys- 
tematically larger than Ay (foreground) for YSOs and possi- 
ble YSOs. Such over all behavior is expected for YSOs, since 
Ay (foreground) from lSchultheis et"aD (120091) is a spatially av- 
eraged line-of-sight extinction to the GC, while Ay (SL+LL) 
is the sum of the line-of-sight extinction to the GC and the 
localized extinction from the dusty envelope of the YSO 

3.4. Molecular Abundances 

By using the dust extinction values derived in the previous 
section, we derived abundances for gas-phase molecular ab- 
sorbers with respect to hydrogen. We obtained a total hydro- 
gen column density from the optical d epth of the 9.7 /nm sili- 
cate absorption, assuming Ay / T9.7 = 9 (Roche & Aitke nl 19*851) 
and N H /A V « 1.87 X 10 21 cm" 2 mag" 1 dBohlin et al.lfl978l) at 
Ry = 3.1. We used Ay(SL+LL) to derive the H2 column den- 
sity, assuming A/h, = Nh/2. Here we implicitly assumed that 
the H2 column density along the full 8 kpc line of sight is 
comparable to the local value near the YSO. A factor of two 
difference would exist, if the local and the full H2 column 
densities are the same, but we neglected this difference. 

The gas-phase molecular abundances relative to H2 (i.e., ra- 
tios of column densities) are shown in Table [5] Our derived 
abundances for C2H2 and HCN are 10~ 6 9 -10~ 5 3 , and our gas- 
phase CO2 abundances are 10~ 6 4 — 10~ 51 . Intervening molecu- 
lar clouds in the line of sight to the GC are less likely the main 
cause of these absorptions because the average HCN abun- 
dance of 2.5 x 10" 8 towards Sgr B2(M) dGreaves & Nvmanl 



1996), where half of our YSOs and possible YSOs are found 
(§ 14.3b . is an order of magnitude lower than our measure- 
ments. 

Individual gas-phase abundances are comparable 
to or generally high e r tha n those in earlier studies. 
Lahui s & van Dishoeckl (2000) found abundances of 10~ 8 - 
10~ 6 for C2H2 and HCN in the warm gas for several massive 
YSOs, and iKnez et all ( f2009h found 10" 6 1 for C 2 H 2 and 
IP" 8 3 for HCN towards IRS 1 in NGC 7538. iBoonman et all 
( 120031) estimated C0 2 abundances of 10" 7 2 -10" 6 5 towards 
lines of sight to several YSOs. However, these differences 
could be due to the uncertainties of comparing different 
techniques of deriving N(H2). If we consider the column 
densities of warm gas towards massive YSOs, our values 
are in good agreement with those found in the previous 
work ( Lahuis & van Dishoeck 2000J IBoonman et all ;2003; 
IKnez et alj|2009l) . In addition, our gas to solid abundance ra- 
tios for CO2 (10 _I -10~ 2 ), which do not requi re knowledge of 
the for eground extinction, are consistent with IBoonman et all 
(12001 . 

4. PROPERTIES OF MASSIVE YSOS IN THE GC 

In the above section, we spectroscopically identified 16 
YSOs and 19 possible YSOs from among 107 IRS targets 
in the GC. Although our selection of massive YSOs is pri- 
marily based on the 15 /im CO2 ice absorption profile, ab- 
sorption from hot and dense molecular gases further supports 
our selection procedures. In this section, we derive and in- 
spect properties of these YSOs and possible YSOs using SED 
model fits, and look for a spatial correlation of these sources 
in the CMZ. 

4. 1 . YSO Parameters from SED Fitting 
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FIG. 13. — Comparisons of (SL +LL) with Av ( SH+LH) (top) and with 
those based on the colors of giant stars in the field (Schultheis et al. 2009) 
(bottom). Red circles are YSOs, and blue triangles are possible YSOs se- 
lected in this work. Dotted line represents equal values for Ay (SL+LL) and 
Ay (SH+LH). Solid line in the top panel shows a linear fit to the data using er- 
rors in both axes. Ay (SL+LL) is systematically larger than Ay (foreground) 
as Ay (SL+LL) is the sum of the line-of-sight extinction and the localized 
extinction from the dusty envelope of a YSO while Ay (foreground) is a spa- 
tially averaged line-of-sight extinction. 



For our 35 YSOs and possible YSOs, we performe d 
SED fitting using a set of models in Robit aille et alj d2006 |). 
For this purpose, we used the Online SED FitterFI 
(Robitaille et al. 2007) to derive YSO parameters, such as the 
mass of the central object, the bolometric luminosity, and the 
accretion rate from the envelope. 

As an input to the SED Fitter, we used available near- 
and mid-IR photometry as listed in Table[7] The near-IR JHK 
observations are Aperture 3 magnitudes from UKIDSS 
DR2 dWarren et all 120071) . Many of our YSOs and possi- 
ble YS Os are found on sa turated pixels on the MIPS 24 /im 
images dCarev et al.l 120091) . Therefore, we derived synthetic 
photometry at 24 by convolving the MIPS [24] filter re- 
sponse function on the IRS spectra, following the prescrip- 
tions on the Spitzer websit43 These values are listed in Ta- 
ble|7] Synthetic values for our IRS targets are 1.05±0.17mag 
(A'comp = 29) systematically smaller (brighter) than MIPS 
[24] photometry (S. Carey, 2008, private communication). 



n |http : //caravan . astro . wise ■ edu/protostars| 

12 See |http : / / ssc . spitzer . caltech . edu/dataanalysistool s/co<51c&6*6K l 



Similarly, we found a mean difference of 0. 72 ± 0.09 mag 
(jVconro = 77) between IRAC [8.0] photometry (iRamfrez et al.l 
2008) and synthetic values. Again the sense of the differ- 
ence is that synthetic values are brighter than Ramirez et 
al. values. This is likely due to extended emission around 
YSOs. We also utilized 450 /im and 850 /im observa- 
tions from th e Submilhmetre Commo n User Bolometer Ar- 
ray (SCUBA; |Pi Francesco et al.ll2008l) . measured with a 23" 
diameter beam. 

In addition to the above photometry, we derived monochro- 
matic fluxes at 14 wavelength points: 5.58 /im, 6.4 /im, 7.65 
/im, 8.5 /im, 9.0 /im, 9.7 /im, 11.0 /im, 12.0 /im, 13.5 /im, 
17.0 /im, 18.0 /im, 21.0 /im, 30.0 /um, and 35.0 /im. These 
points were selected to characterize the overall shape of a 
SED with as little ice features as possible, because the models 
do not include ices. We computed a monochromatic flux with 
a 2%-wide Gaussian filter in these wavelength points, except 
at 9.7 /im where we used a 3%-wide filter, to better character- 
ize the bottom of the silicate absorption band. Note that we 
did not use the IRAC [5.8] and [8.0] photometry and instead 
used the above synthetic values to avoid strong 6 /im and 7 /im 
absorption bands, which are not included in the models. 

For each source, we ran the Online SED Fitter us- 
ing the above set of photometry and collected results that sat- 
isfy (x 2 -xLi)Mot < 5, where x^ n is the minimum x 2 value 
from the available model sets, and N tot is the total number of 
data points, which are between 12 and 22 for our sources. We 
note that the fitting is not strictly statistical, given the lim- 
ited p^ame^e£_sgaceofmodels for all 14 YSO parameters 
(see Robita ille et al.120 07). We chose the above cut to include 
reasonable fitting results, and then estimated a mean and a 
standard deviation for each YSO parameter. 

Figure [14] displays SED fitting results overplotted on the 
input photometry for four YSOs. All 35 YSOs and possible 
YSOs are shown in the online journal. The black solid line 
shows a best-fitting SED, and grey lines show acceptable fits 
from the on-line SED fitter. The dashed line represents the 
emission from the central source in the absence of extinction 
from the dusty envelope. In this fitting exercise, we restricted 
the source distance, d, to 7 kpc < d < 9 kpc from the Sun, 
and interstellar extinction along the line of sight to the GC to 
20 mag < Ay < 40 mag. 

Figure [15] shows results for derived YSO parameters, and 
Table [8] summarizes the results. Entries with no error bars in- 
dicate that a single solution is found within (x 2_ Xmin)/^tot < 
5. Our derived masses of central objects span 8 Mq < M* < 
23 Mq, and the total luminosities range over 10 3 L© < L tot < 
10 5 L©. Note that the mass is not directly determined from 
the SED; rather, it is the bolometric luminosity and the tem- 
perature we are determining, and the mass is implicitly con- 
strained by these from the evolutionary tracks built-in to the 
model grid. The mass accretion rate from the envelope is 
another indicator for the evolutionary stage of YSOs. For 
our YSOs and possible YSOs, we found a heavy infall rate, 
10" 4 M Q yr~ l < M env < 1 0~ 3 M(p yr" 1 , which is consistent with 
those for Stage-I YSOs dRobitaille et al.ll2006l) . The range of 
these parameters recovered from the SED fitting tool 
remained essentially unchanged if we instead imposed a Ay 
limit using Ay (foregro und) measurements i n Table |6] with 
its ±2a error bounds (iSchulfheis et al.l 12009b . Our SED fit- 
ting suggests that our sources are massive YSOs in their early 
stages of protostar evolution. 



i ^l s Vi 0WS the color distribution of YSOs and possible 
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FIG. 14. — SED fitting results for four GC YSOs using a set of models in Robitaille et al. (2006) and based on near-infrared photometry, synthetic values 
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the dusty envelope. The SED model fitting suggests that our spectroscopically selected YSOs are massive Stage-I YSOs. (An extended version of this figure is 
available in the online journal, showing all 35 GC YSOs and possible YSOs.) 



YSOs in the mid-IR color-color diagrams, o verlaid with re- 
gions occupied by theoretical Stage-I objects (Robitaille et al. 
2006). The colors of YSOs and possible YSOs, relative to 
non-YSOs, are a bit bluer for [5.8]-[8.0]. On the other hand, 
YSOs and possible YSOs are redder in [3.6]-[4.5] and [3.6]- 
[5.8]. All of our 107 sources have similar [8.0]-[24] colors. 
Although YSOs and possible YSOs in the GC have colors that 
are similar to the theoretically predicted colors, non-YSOs 
are also found in the s ame color space. Thi s confirms earlier 
theoretical work (e.g., Robitaille et al. 2007), concluding that 
broad-band colors are not sufficient to separate YSOs from 
non-YSOs. 

4.2. Mass Estimates from Radio Continuum 

Eight YSOs and possible YSOs are coincident with ra- 
dio continuum sources, and are thus likely to be com- 
pact H II regions. These are listed in the last two 
columns of Table [H We used radio continuum data 
dZoonematkermani et all 1 19901; iMehringer et all 1 19921 1 1993 ; 
Becker et al.lll994t iMehringerl 119951; lLazio & Cordesl Il998[ 
Yusef-Zadeh et dmoOTlWhite et all 120051: lLazio & Cordesl 
2008; Yusef-Zadeh et al. 2009) to derive the number of ioniz- 
ing photons for each H II region, assuming a distance of 8 kpc. 
We then converted the numbe r of ion i zing p hotons to stellar 
mass by using the results of Panagia (1973) and assuming a 
surface gravity of logg = 4.2. Our derived masses, listed in 
Table [8j agree on average with those estimated from SED fits 
in the previous section (§ 14.1 I t. 

4.3. Spatial Distribution of YSOs in the GC 



We confined our spectroscopic sample to those within \b\ < 
15' to avoid likely foreground objects (§ 12. It . Nevertheless, 
this spatial cut is generous enough that our spectroscopic sur- 
vey is almost free from a spatially dependent sample bias, and 
enables us to map out active star-forming regions in the GC 
and to study their relation to the interstellar medium (ISM). 

Figure [T71 displays the locations of 35 YSOs and possible 
YSOs in the CMZ (see Figure [TJfor the locations of all of our 
spectroscopic targets). Although YSOs and possible YSOs 
are found throughout the CMZ, it is striking to see that half 
of these sources (18 out of 35) are found in and around Sgr B. 
Sgr B is known as the mo st active star-forming region in the 
Galaxy dBallv et alj |2010). but this is the first direct evidence 
of the presence of YSOs in this region at the earliest stage of 
star formation (< 1 Myr). Figure [18] shows the Sgr B region 
with locations of our YSOs and possible YSOs. As seen on 
the 24 ixm map, our sources are preferentially found on the 
edge of strong 24 /im emission regions. 

4.4. Star Formation Rate at the GC 

YSOs are direct tracers of early star formation, and can be 
used to estimate the in situ star formation rate (SFR) in the 
GC. Previous identifications of YSOs based on broad-band 
photometry were used to infer the SFR in the GC, but the 
heavy extinction towards the GC limits any estimate of the 
SFR based on photometrically selected YSOs. This is because 
reddened AGB stars have similar colors (e.g., Figure [TBI. Our 
IRS spectra provide a unique opportunity to check how well 
earlier studies selected their YSO candidates, and can be used 
to refine SFR estimates at the GC. 
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Table [2] includes cross-identifications of our I RS sample 
with e a rlier photometry-base d YS O selections in IFelli et al] 
(2002). ISchuller et all d2006h. andlYusef-Zadeh et ^"( 2009). 
YSO selections in both IFelli et all (120021) and ISchuller et alJ 
(120061) are based on ISOGAL phot ometry (lOmont et al.|2003b 
ISchuller et ai]|2003l) . while that o f lYusef-Zadeh et al.1 ho6% 
is based on the Spitzer IRAC dRamfrez et all 120081) and 
MIPS photometry (JHinz et al. 20 Q9|). The source catal ogs 
(lOmont et alj|2003t ISchuller et alj |2003: Hin z et al.1120091) do 
not cover the entire CMZ, in particular near Sgr A. Sources 
with missing data (" ••• ") in Table [2] represent our spectro- 
scopic targets that were not detected in these catalogs in a 3" 
s earch radius 

IFelli etal J (120021) used ISOGAL photometry at 7 fim and 15 
fim to select bright YSO candidates, using the mid-infrared 
color-magnitude diagram for ultra-compact H II regions. In 
total, 28 sources identified by IFelli et al.l (120021) as photomet- 
ric YSOs ("yes" in column 6 of Table |2|i were cross-matched 
with our IRS targets (Table [2] column 5) in a 3" search ra- 
dius, but we identified only 36% of them (10/28) as YSOs (4) 
or possible YSOs (6) in our study. 

YSO candidates were also selected by Sch uller et al.1 (120061) 
based on ISOGAL photometry at 7 /im and 15 fim and spa- 
tial extent of ISOGAL sources. Their study focused on a 
small 20' x 20' field between Sgr A and Sgr C. We have 
obtain ed IRS spectra of o nly eight ISOGAL sources in this 
field. ISchuller et ail (120061) photometrically identified five GC 
sources as YSOs (Tabled column 7). However, none of them 
are identified by us as YSOs or possible YSOs (Table [2] col- 
umn 5). The low rate in the YSO identification could be 
due to their selection criteria based on the spatial extent of 



sources, while our spectroscopic targets wer e selected from 
point sources in the IRAC bandpasses (§ 12. 1|. 

A comparable hit rate to that from Felli et al. (20021) was 
found for YSO candidates from the most recent photometric 
study by lYusef-Zadeh et alJ d2009l) . whose YSO candidates 
were identified based on the Spitzer IRAC and MIPS images. 
In total 17 photometric YSOs ("yes" in Table [2] column 8) 
in their list were cross-matched with our IRS targets in a 3" 
search radius, but only 47% (8/17; "yes" in Table [2] column 
8) of them were found to be either YSOs (3) or possible YSOs 
(5) in our study. 

A complete analysis on the SFR estimate requires a bet- 
ter understanding of the sample bias in our spectroscopic tar- 
get selection, which is the subject of the next papers of this 
series. Nonetheless, we can make a preliminary estimate 
on the SFR based on the result in this p aper: since the hit 
rate of the photometric YSO selection in lYusef-Zadeh et al.1 
(2009) is ~ 50%, their SFR estimate for Stage I YSOs would 
have been overestimated by a factor of ~ 2. They have con- 



cluded that the Stage I SFR is 



0.14 M Q yr so this im- 



plies a revised SFR ~ 0.07 M yr' 1 at the GC. If we as- 
sume a gas su rface density of the GC from the total mass of 
5.3 x iff M Q dPierce-Price et al.ll2000h over the entire CMZ, 
both values of the star formatio n rate are roughl y consistent 
with the Kennicutt-Schmidt law jKennicutI T998). 



5. SUMMARY 

We obtained Spitzer/IRS spectra for 107 sources in the GC, 
which were selected based on near- and mid-IR photometry 
including those obtained from Spitzer/IRAC. Based on the 
shape of the 15 /im CO2 spectral feature and the strength 
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of the 15.4 fim shoulder CO2 ice component, we selected 35 
YSOs and possible YSOs. Our identifications are further sup- 
ported by the presence of hot and dense gas-phase molecular 
absorptions such as C2H2, HCN, and CO2 for some YSOs. 
This is the first spectroscopic identification of a large YSO 
population, tracing an early stage of star formation in the GC. 
Spectroscopic confirmation of candidate YSOs in the GC is 
essential because the older stellar population in the GC, when 
reddened by Ay ~ 30, has infrared colors similar to those of 
YSOs. 

From the SED model fitting, we inferred that the masses of 
these objects are typically ^8-23 M Q , and that the high in- 
fall rate from the envelope suggests that they are on Stage I, 
an ear ly evolutionary stage of protostars (e.g., Robitaille et al. 
120061) . We found that these YSOs and possible YSOs are 
found throughout the whole CMZ, but half of them are located 
in and around the Sgr B. We found that about 50% of photo- 
metrically selected YSOs are spectroscopically confirmed by 
our study. We estimated a preliminary s tar formation rate, 
based on an earlier photometric study by lYusef-Z adeh et al. 
(2009), to be - 0.07 M Q yr" 1 . 

Our SpitzerflRS survey is limited to YSOs of at least 
~ 3 Mq (masses of central objects). However, next gen- 
eration telescopes, such as the Giant Magellan Telescope 
(GMT) or the James Webb Space Telescope (JWST), will 
overcome this limit, exploring significantly less massive stars 
with high-resolution imaging and moderate/high-resolution 
spectroscopic capabilities in the near- and mid-IR range, al- 
lowing detailed studies of the initial mass function in these 
crowded fields. Until then, our SpitzerflRS data will remain 
as a unique database for studying the star formation process 
in the GC. 
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TABLE 1 
IRS Sample of Candidate YSOs 



Source ID 
SSTGC 


R.A. 
(J2000.0) 


Decl. 
(J2000.0) 


Date of 
Observation 






YSO Status 




This work" 1 


H'elli et al. 


Schuller et al. 


Yuset-Zadeh et al. 


244532 


17 43 47.97 


-29 38 41.2 


Oct. 2008 


no 


no 




no 


260956 


17 43 55.98 


-29 36 22.4 


Oct. 2008 


no 


yes 






263857 


17 43 57.32 


-29 36 40.6 


Oct. 2008 


no 




yes 


284291 


17 44 06.91 


-29 24 17.4 


May 2008 


no 


yes 




yes 


293528 


17 44 11.20 


-29 26 37.9 


May 2008 


no 


yes 




300758 


17 44 14.49 


-29 23 22.2 


May 2008 


maybe 


yes 




yes 


303865 


17 44 15.85 


-29 20 43.7 


May 2008 


no 


no 






304239 


17 44 16.03 


-29 33 16.6 


Oct. 2008 


yes 


no 






343554 


17 44 31.54 


-29 27 39.0 


Oct. 2008 


no 


yes 




yes 


348392 


17 44 33.41 


-29 27 02.0 


Oct. 2008 


no 


yes 




349071 


17 44 33.68 


-29 13 55.7 


May 2008 


no 


yes 


yes 




354683 


17 44 35.87 


-29 27 44.8 


Oct. 2008 


no 








358370 


17 44 37.26 


-29 28 41.7 


Oct. 2008 


no 








360055 


17 44 37.90 


-29 25 46.5 


Oct. 2008 


no 


yes 






360559 


17 44 38.09 


-29 28 38.9 


Oct. 2008 


maybe 






368854 


17 44 41.29 


-29 24 35.4 


May 2008 


no 


no 






370438 


17 44 41.90 


-29 23 32.2 


May 2008 


maybe 


no 




yes 


372630 


17 44 42.79 


-29 23 16.3 


May 2008 


maybe 


no 




yes 


374813 


17 44 43.59 


-29 20 48.8 


May 2008 


no 


no 




yes 


381931 


17 44 46.32 


-29 27 39.3 


Oct. 2008 


no 






yes 


388790 


17 44 48.94 


-29 23 42.8 


Oct. 2008 


no 


no 




394248 


17 44 51.02 


-28 50 46.6 


May 2008 


no 


no 




no 


395805 


17 44 51.68 


-29 11 00.2 


May 2008 


no 








401264 


17 44 53.73 


-29 23 12.5 


Oct. 2008 


no 






yes 


404312 


17 44 54.89 


-29 14 13.1 


Oct. 2008 


no 


yes 


yes 


405235 


17 44 55.25 


-29 15 37.8 


Oct. 2008 


no 


yes 


yes 




412509 


17 44 58.01 


-29 10 56.6 


May 2008 


no 




no 


yes 


421092 


17 45 01.27 


-29 14 55.7 


Oct. 2008 


no 




no 


yes 


425399 


17 45 02.91 


-29 22 11.2 


Oct. 2008 


no 


no 




no 


426214 


17 45 03.21 


-29 17 38.3 


Oct. 2008 


no 


yes 


yes 




440424 


17 45 08.58 


-28 46 17.7 


May 2008 


no 


yes 


no 


465659 


17 45 18.10 


-29 04 40.6 


Oct. 2008 


no 


yes 


yes 




492222 


17 45 27.95 


-28 56 22.7 


May 2008 


no 




496149 


17 45 29.42 


-29 10 21.8 


Oct. 2008 


maybe 


no 


no 


yes 


497500 


17 45 29.91 


-28 54 22.8 


May 2008 


no 






507261 


17 45 33.50 


-28 54 37.2 


May 2008 


no 








511261 


17 45 34.94 


-29 25 10.3 


Oct. 2008 


no 








511666 


17 45 35.08 


-28 53 34.2 


May 2008 


no 








516435 


17 45 36.84 


-28 52 21.2 


May 2008 


no 








516756 


17 45 36.94 


-28 54 33.4 


Oct. 2008 


no 








517724 


17 45 37.30 


-28 53 53.7 


May 2008 


no 








519103 


17 45 37.80 


-28 57 16.2 


May 2008 


no 








521894 


17 45 38.82 


-28 52 31.9 


May 2008 


no 








524419 


17 45 39.80 


-28 53 44.4 


May 2008 


no 








524665 


17 45 39.86 


-29 23 23.4 


Oct. 2008 


yes 








525666 


17 45 40.22 


-28 53 28.2 


May 2008 


no 









a YSOs are marked as "yes", possible YSOs are marked as "maybe", and the remaining targets are marked as "no". 
b 4.5 fim excess source without a 24 (im counterpart. 
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TABLE 2 
IRS Sample of Candidate YSOs 



Source ID 


R.A. 


Decl. 


Date of 




YSO Status 


SSTGC 


(J2000.0) 


(J2000.0) 


Observation 


This work' 1 


helli et al. Schuller et al. Yusel-Zadeh et al. 


531300 


n ac ia oo 

17 45 42.32 


O O CO All 

—28 52 47.3 


a it onno 

May 2008 


no 




534806 


n A c Al C"7 

1/45 43.57 


o o on i c n 

—28 29 16.9 


/~i . . onno 

(Jet. 2008 


no 


no 


53500/ 


n AC Al 

1/45 43.64 


o o co o a n 

—28 52 24.9 


a x onno 

May 2008 


no 




536969 


17 45 44.35 


o n / vi 1 o o 

-29 01 13.8 


/~\ i. OftftO 

Oct. 2008 


no 




C A f\ O A n 

540840 


n ac i a 

1/45 45. /4 


oo a o on i 

28 48 29. / 


a x onno 

May 2008 


no 




C A ICCi 1 

543691 


1/45 46. /6 


on no /ion 
—29 02 48.0 


c\,^ onno 
UCt. 2008 


no 




C A 7 O 1 7 

54/81 / 


n if aq o /i 
1/ 45 48.24 


O O A O 1 C C 

—28 48 16.6 


onno 
May 2008 


no 




550608 


n a c /in in 
1/ 45 49.30 


oo cn coo 
—28 50 58.8 


c\^t onno 
UCt. 2008 


no 




56 J 780 


in AC CA 11 

1/45 54.1 1 


O O CO 1 O 1 

—28 58 12.1 


/-\ , onno 

(Jet. 2008 


maybe 




564417 


n C/1 ii 

1/45 54.33 


on nn t\i o 
—29 00 03.2 


/-\ , onno 
(Jet. 2008 


no 




579667 


1 1 ac c n nf\ 

17 45 59.90 


o o Co f\i o 

—28 53 07.2 


Oct. 2008 


no 




f Ofl1 O O 

580183 


1 —J AC (\f\ f\1 

1/46 00.0 / 


—29 01 49.3 


/■ v . , onno 

(Jet. 2008 


no 


yes 


584613 


n A C A1 ci 

LI 46 01.6 / 


OO QC ci n 

—28 35 53.9 


"Nifn.r onno 
May 2008 


no 


no 


588220 


11 ac no no 

17 46 02.98 


o o co a c n 

—28 52 45.0 


/ v . , onno 

(Jet. 2008 


no 




6002 / 4 


1 i a c m q n 
1/46 o/.3y 


o o a c io n 
—28 45 32.0 


»(„,, onno 
May 2008 


no 




609613 


1/46 10. / 1 


oo /i o cc n 
—28 48 55.0 


onno 
(Jet. 2008 


no 




6 1 0642 


n ;i 11 no 
1 / 46 1 1.08 


oo cc /in n 
—28 55 40.9 


»«„,, onno 
May 2008 


maybe 




612688 


1 T A £. 1 1 Ol 

1/46 11.83 


O O A 1 11 

—28 4/ 12.0 


May 2008 


no 




c 1 on 1 o 
6 1 80 1 8 


1/46 13.81 


OO A 1 A A C 

—28 43 44.5 


ni,„ onno 
May 2008 


maybe 


yes 


619522 


n /i c i a ii 
1/46 14.33 


O O A 1 1 O /I 

—28 43 18.4 


»«„,, onno 
May 2008 


maybe 


C 1 CittCA 
619964 


1/46 14.48 


oo i/^ an 7 
—28 36 39. / 


»(„,, onno 
May 2008 


no 


no ■ ■ ■ no 


£7 lOCO 

621858 


1/46 15.18 


OO CO 11 A 

—lis 52 31.4 


/-»„f onno 
(Jet. 2008 


no 




635358 


n a c on n i 

17 46 20.01 


O O ACt 1 O 1 

—28 49 18.3 


/-\ , onno 

(Jet. 2008 


no 




64602 1 


17 /l A tj on 
1/ 46 23. 8y 


70 1Q A Q 1 
—28 3y 48. 1 


A,fn7 onno 
May 2008 


no 


no 


AA 07On 

648 /yo 


1 H AC. 7/1 oi 
1 / 46 24. y3 


70 A 7 10 O 
—28 4 / 18.2 


o^t onno 
UCt. 2008 


no 




6532/0 


1/ 46 26.55 


o o 1 o cn n 

—28 18 5y.y 


r~\ „f onno 
UCt. 2008 


maybe 


yes • • • no 


660/08 


it /I/: on 07 
1 / 46 2y.2 / 


OQ CA (\1 O 

—28 54 03. y 


a^h7 onno 
May 2008 


no 


no ■ ■ ■ 


*nnnco 

6/0953 


11 AC TO (\ C 

17 46 32.95 


o o a o i /; i 

—28 42 16.3 


May 2008 


maybe 


no 


673 15 1 


17 AC 11 1 C 

1/ 46 33. /6 


o o /in i o n 
—28 40 32. y 


»«„,, onno 
May 2008 


no 


no ■ ■ ■ no 


C1C\C\1C 

6/y036 


\i ac ic no 
1/ 46 35. y8 


OO Al CO O 

—28 43 58.2 


»(„,, onno 
May 2008 


maybe 


yes 


68y3y / 


mac in CI 
1/ 46 3y.6 / 


OO A 1 07 O 

—28 41 2 /.8 


ni,„ onno 
May 2008 


no 


yes 


Ct\C1C1 

6y636/ 


17 AC /lO OO 

1/46 42.28 


O O Q 1 HC 1 

—28 33 26.3 


»«„,, onno 
May 2008 


no 


no ■ ■ ■ no 


71 1462 


1 1 AC Al OO 

1/46 4/. 82 


OO Al 1 C A 

—28 4 / 15.4 


»«„,, onno 
May 2008 


no 


yes • • • no 


716531 


11 AC /in CA 

1/46 49.64 


OO 1C CI A 

—28 36 5 /.4 


/-\ , onno 

(Jet. 2008 


no 


yes 


1 1 C7C7 

/18/5 / 


ii ac cn cn 

17 46 50.50 


OO /II 11 A 

—28 43 33.4 


* jt onno 

May 2008 


maybe 


yes 


719445 


11 A /T c r\ TO 

17 46 50.72 


OO 11 i~\ A 1 

—28 31 24.7 


May 2008 


yes 


yes • • • yes 


721436 


n ac. ci ac\ 
1/46 51.49 


lO 11 (\C o 

—28 33 06.2 


a jt . onno 

May 2008 


no 


no 


722141 


11 AC CI CQ 

1/46 51.68 


OO OO A 1 C 

28 28 41.6 


a it onno 

May 2008 


yes 


yes • • • yes 


/2632 / 


11 ac ci on 
1/ 46 53. 2y 


o o i o n 1 o 
—28 32 01.2 


c\„ t onno 
(Jet. 2008 


yes 


yes 


/28480 


1 1 AC CA A 1 

1/ 46 54. 13 


oo on an c 
—28 2y 3y.5 


ni,„ onno 
May 2008 


yes 


yes 


73253 1 


11 AC ?C 1 A 

1/ 46 55. /4 


oo qo on o 
—28 32 20.2 


r~\ „f onno 
(Jet. 2008 


no 




•70 on/: 

/38126 


1 1 AC ZI C\C 

1 / 46 5 / .y5 


OO IC C/1 c 

—28 35 54.5 


»«„,, onno 
May 2008 


no 




/ 606/9 


11 Al (\1 A C 

17 47 07.45 


o o oo a 1 n 
—28 28 41.9 


a t . onno 

May 2008 


yes 




761771 


ii ah m n/i 

1 / 4 / o/.y4 


OO O/l CQ O 

—28 24 53.2 


ni,„ onno 
May 2008 


yes 




/oyjuj 


1 "7 /1 7 11 97 
1/4/ 11. LI 


70 7/^ 1~\ 7 
— ZO ZD J 1 . / 


A/Tni7 7nno 
iviay zuuo 


yes 


no ■ ■ ■ 


770393 


17 47 11 75 


—98 31 91 Q 

ZO JI Z 1 . 7 


fVf 9008 


yes 




771791 


17 47 1 ? 35 


—98 31 10 8 


Oct 9008 


no 




772151 


1 7 47 1 ? 50 


—98 94 1 5 6 

Z O Z'T 1J.U 


Mav 9008 


yes 




77900 i 


1 7 47 19 QO 


—98 39 05 5 
ZO jZ Uj.j 


Clrt 9008 


yes 




773985 


17 47 13.34 


-28 31 56.9 


Oct. 2008 


maybe 




782872 


17 47 17.31 


-28 32 20.2 


Oct. 2008 


no 


yes 


786009 


17 47 18.69 


-28 27 31.7 


May 2008 


maybe 


790317 


17 47 20.55 


-28 23 54.8 


May 2008 


maybe 




797384 


17 47 23.68 


-28 23 34.6 


May 2008 


yes 




799887 


17 47 24.80 


-28 15 56.8 


May 2008 


maybe 


yes • • • no 


801865 


17 47 25.69 


-28 24 40.2 


May 2008 


yes 




803187 


17 47 26.29 


-28 22 01.5 


May 2008 


yes 


yes b 


803471 


17 47 26.40 


-28 24 43.7 


May 2008 


yes 




806191 


17 47 27.66 


-28 26 28.4 


Oct. 2008 


maybe 


no • • • yes 


817031 


17 47 32.97 


-28 34 12.0 


Oct. 2008 


no 


yes • • • no 


817663 


17 47 33.28 


-28 24 47.4 


May 2008 


maybe 



a YSOs are marked as "yes", possible YSOs are marked as "maybe", and the remaining targets are marked as "no". 
b 4.5 fim excess source without a 24 fim counterpart. 
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TABLE 3 

Data Reduction Summary for YSOs and possible YSOs 



SSTGC 


Baseline 




Flux Scaling Factor b 




Excluded 


ID 


Module 11 


SL1 


SL2 


SL3 


LLI 


LL2 


Data Set 


1 A A "7 C O 

300/35 


LL2 


1 OC _L_ A AO 

1.25 ± 0.02 


1.18±0.06 


1.15 ±0.01 


1.03 ±0.03 






304239 


SL1 




0.89 ±0.08 


1.31 ±0.02 


1.80±0.14 


1 O 1 I A 1 A 

l.8l ±0.14 


SL2 (1st nod, N), LL (W) 


360559 


SL1 




2.34 ±0.06 


2.70 ±0.14 


0.10±0.01 


A 1 A 1 A A 1 
0.10 ±0.01 




3 /043o 


LL2 


1 in _i_ a 1 n 
1.30 ± 0. 10 


1.17 ±0.08 


1.23 ±0.03 


0.98 ±0.04 






3 / 2630 


LL2 


1 Q C _L a AT 

1.85 ± 0.0/ 


1.60 ±0.15 


1.61 ±0.02 


0.96 ±0.03 




5L1 (£.), bLZ (2nd nod, JN, r.) 


A C\£L 1 A A 

496149 


LL2 


("1 0"7 _L_ ("1 AT 

0.87 ± 0.03 


0.85 ±0.04 


0.92 ±0.01 


1.00 ±0.02 






CO/1 

524665 


SL1 




1.03 ±0.06 


1.06 ±0.01 


0.86 ±0.01 


A OO _1_ A AO 

0.52 ± 0.02 




563 /80 


SL1 




1.04 ±0.05 


1.06 ±0.01 


0.27 ±0.01 


A O C _l_ A A 1 

0.25 ± 0.01 


SL1 (W) 


610642 


LL2 








0.97 ±0.02 




SL not available 


618018 


LL2 


1 1 A AO 

1.32 ±0.02 


1.34 ±0.06 


1.34 ±0.01 


1.02 ±0.02 






619522 


SL1 




1.05 ±0.03 


1.03 ±0.02 


0.54 ±0.02 


A C O _l_ A AO 

0.58 ± 0.02 




ACQOOA 

6332 /0 


LL2 


0.99 ± 0.01 


0.99 ±0.07 


0.99 ±0.01 


1.02 ±0.03 






"7 a a c ^ 

6/0953 


LL2 


O A£ _l_ A 1 1 

2.06 ± 0.1 1 


2.09 ±0.06 


2.10±0.02 


1.00 ±0.03 




SL2 (1st nod) 


6 /9036 


LL2 


1 1£ _L A 1 A 

1.36 ± 0. 10 


1.33 ±0.08 


1.36±0.01 


0.94 ±0.04 




SL1 (W) 


O" 1 0"7C"7 


LL2 


1 _l_ A AO 

1.66 ± 0.02 


1.56 ±0.09 


1.57 ±0.02 


0.98 ±0.02 




SL2 (E) 


1 1 9445 


LL2 


1 on I A A A 

1.29 ± 0.04 


1.25 ±0.09 


1.28 ±0.01 


0.97 ±0.02 






722141 


LL2 


1.81 i 0.1 1 


1.77±0.11 


1.85 ±0.02 


0.96 ±0.02 




SL1 (W) 


/2632 / 


LL2 


1.81 ±0.11 


1.77±0.11 


1.85 ±0.02 






bLl (W), bLZ (1st nod, JN), LL1 (both nods) 


TO O /I OCA 

/28480 


LL2 


O 1 1 I A A A 

1.1 \ ± 0.04 


2.60 ±0.18 


2.69 ±0.02 


1.01 ±0.02 




SL1 (W) 


/ 606 79 


LL2 


1 CI _L_ A A A 

1.53 ± 0.04 


1.50 ±0.05 


1.51 ±0.02 


0.95 ±0.01 




SL1 (E) 


761771 


LL2 


1 1 _L_ A CA 

3.21 ± 0.50 


2.93±0.16 


3.25±0.13 


0.98 ±0.02 






/69305 


LL2 


1 A O _L A AA 

1.42 ± 0.06 


1.39 ±0.09 


1.42 ±0.02 


1.00 ±0.04 






/ /03V3 


LL2 


-1 OO _l_ A A/1 

3./y ± 0.04 


3.12±0.18 


3.28 ±0.02 






lli (Dotn noasj 


772151 


LL2 


1.37 ±0.08 


1.35 ±0.05 


1.54 ±0.03 


1.01 ±0.06 




SL1 (S), SL2 (W) 


772981 


LL2 


2.20 ±0.08 


1.88±0.01 


1.99 ±0.01 


0.64 ±0.02 




SL2 (1st nod, N), LLI (N) 


773985 


LL2 


7.87 ±1.59 


6.02±0.10 


6.51 ±0.14 


0.90 ±0.03 




LLI (N, S) 


786009 


LL2 


1.45 ±0.03 


1.81 ±0.13 


1.89 ±0.06 


1.00 ±0.02 




SL1 (1st nod, S) 


790317 


LL2 


1.45 ±0.05 


1.45 ±0.09 


1.49 ±0.01 


0.93 ±0.02 




SH(S) 


797384 


LL2 


1.33 ±0.05 


1.34 ±0.09 


1.39 ±0.01 


0.99 ±0.03 






799887 


LL2 


1.02 ±0.06 


1.02 ±0.07 


1.00 ±0.01 


0.96 ±0.02 






801865 


SL1 




0.93 ±0.07 


1.00 ±0.01 


0.05 ±0.01 


0.04 ±0.01 




803187 


LL2 


1.33 ±0.13 


1.31 ±0.08 


1.33 ±0.01 


0.98 ±0.03 






803471 


LL2 


1.83 ±0.04 


1.85±0.13 


1.90±0.01 


1.03 ±0.03 






806191 


LL2 


1.37±0.14 






1.10 ±0.03 




N from all sky positions, SL2 (W), SL3 (both nods, S) 


817663 


SL1 




0.69 ±0.04 


0.72 ±0.05 


0.42 ±0.01 


0.48 ±0.01 


SL2 (E) 



Note. — IRS modules: short-high (SH; 9.9 ^im-19.6 (im, A/AA ~ 600), long-high (LH; 18.7 /im-37.2 (im, A/AA ~ 600), short-low (SL [1st order SL1 7.4 ^im-14.5 fj,m, 2nd 
order SL2 5.2 (im-7.7 fim, 3rd order SL3 7.3 (im-8.7 (im, A/AA ~ 60-127]), and long-low (LL [1st order LLI 19.5 ^m-38.0 fim, 2nd order LL2 14.0 (im-21.3 (im, 3rd order 
LL3 19.4 ^im-21.7 (im, A/AA ~ 57-126]). 



a IRS module selected as a baseline for the flux calibration. See text. 

k Adopted scaling factor in each module. The value represents the scaled flux divided by the original flux. 

c Specific modules/orders that contain defective data. These were excluded in the spectral analysis. "NSEW" denote background observations northern/southern/eastern/western from 
the source target. "Nod" represents a specific nod position for a set of target spectra. 
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TABLE 4 

C0 2 Ice Decomposition for YSOs and possible YSOs 



SSTGC 


log N col (polar) 


logW C oi(apolar) 


log W co i (shoulder) 


log A^oi (diluted) 


log N co] (pure) 


log Wei (total) 


Xtot 


Mot 


YSO 


ID 


(cm- 2 ) 


(cm- 2 ) 


(cm- 2 ) 


(cm- 2 ) 


(cm" 2 ) 


(cm- 2 ) 






Status 


300758 


17.08 ±0.08 


17.11 ±0.06 


16.88 ±0.06 


< 16.17 


< 16.78 


17.57 ±0.03 


488 


152 


maybe 


304239 


18.31 ±0.08 


17.98 ±0.06 


17.40±0.16 


< 16.17 


17.29 ±0.25 


18.54 ±0.04 


1615 


152 


yes 


360559 


17.40±0.17 


17.41 ±0.15 


16.88±0.21 


16.54±0.31 


< 16.78 


17.79 ±0.01 


2181 


152 


maybe 


370438 


17.83 ±0.19 


17.42 ±0.05 


17.16±0.08 


< 16.17 


< 16.78 


18.06±0.12 


280 


152 


maybe 


372630 


17.06 ±0.16 


17.48 ±0.01 


17.11 ±0.05 


< 16.17 


< 16.78 


17.73 ±0.06 


237 


152 


maybe 


496149 


17.70±0.10 


17.42 ±0.04 


16.97 ±0.11 


< 16.17 


< 16.78 


17.95 ±0.04 


99 


152 


maybe 


524665 


17.75 ±0.01 


17.41 ±0.05 


17.16±0.02 


< 16.17 


16.94 ±0.06 


18.02±0.01 


424 


152 


yes 


563780 


17.06±0.01 


17.05 ±0.07 


16.75 ±0.06 


< 16.17 


< 16.78 


17.39±0.11 


594 


152 


maybe 


610642 


17.31 ±0.06 


16.80 ±0.09 


16.79 ±0.04 


< 16.17 


< 16.78 


17.54 ±0.03 


219 


118 


maybe 


618018 


17.18±0.03 


17.10±0.03 


16.86 ±0.03 


< 16.17 


< 16.78 


17.55 ±0.02 


535 


152 


maybe 


619522 


17.60 ±0.69 


17.19 ±0.30 


17.07 ±0.23 


< 16.17 


< 16.78 


17.83 ±0.62 


1027 


152 


maybe 


653270 


17.99 ±0.05 


16.79±0.21 


16.89±0.11 


16.82 ±0.23 


17.39±0.17 


18. 14 ±0.04 


980 


152 


maybe 


670953 


17.54 ±0.39 


16.90 ±0.30 


16.78 ±0.20 


16.64 ±0.31 


< 16.78 


17.75 ±0.14 


710 


152 


maybe 


679036 


17.26 ±0.15 


17.17 ±0.04 


16.90 ±0.06 


< 16.17 


< 16.78 


17.62 ±0.04 


473 


152 


maybe 


718757 


17.76 ±0.03 


17.59 ±0.01 


16.96 ±0.01 


16.52 ±0.04 


< 16.78 


18.04 ±0.02 


658 


152 


maybe 


719445 


17.88 ±0.03 


17.32 ±0.04 


17.14±0.03 


16.36±0.15 


< 16.78 


18.07 ±0.03 


889 


152 


yes 


722141 


17.80 ±0.02 


17.48±0.01 


17.12 ±0.02 


16.24 ±0.12 


< 16.78 


18.03 ±0.01 


990 


152 


yes 


726327 


17.47 ±0.10 


17.30 ±0.07 


17.08 ±0.07 


< 16.17 


< 16.78 


17.80 ±0.01 


867 


152 


yes 


728480 


17.83 ±0.01 


17.32 ±0.01 


17.04 ±0.01 


16.20 ±0.03 


< 16.78 


18.00 ±0.01 


1175 


152 


yes 


760679 


17.70 ±0.01 


17.34 ±0.02 


17.11 ±0.03 


16.23 ±0.06 


< 16.78 


17.94 ±0.01 


834 


152 


yes 


761771 


18.24 ±0.03 


17.73 ±0.01 


17.33 ±0.02 


16.42 ±0.05 


< 16.78 


18.40 ±0.02 


963 


152 


yes 


769305 


18.21 ±0.12 


18.03 ±0.06 


17.57 ±0.05 


16.68 ±0.04 


< 16.78 


18.49 ±0.04 


2347 


152 


yes 


770393 


17.48 ±0.03 


17.23 ±0.01 


16.78 ±0.01 


< 16.17 


< 16.78 


17.74 ±0.02 


1077 


152 


yes 


772151 


18.13 ± 0.10 


1"7 QA _1_ O IIA 

1 /.o4 ± 0.06 


1 "7 A 1 _L n m 
1 /.41 ± 0.0/ 


16.56 ± 0.31 


< 16. /o 


lo.3 / ± O.Uj 


1081 


152 


yes 


772981 


17.15±0.15 


17.19 ±0.07 


16.95 ±0.05 


< 16.17 


< 16.78 


17.60±0.13 


630 


152 


yes 


773985 


17.26 ±0.14 


17.24 ±0.03 


16.85 ±0.05 


< 16.17 


< 16.78 


17.64 ±0.04 


340 


152 


maybe 


786009 


17.43 ±0.02 


17.16 ±0.02 


16.81 ±0.02 


< 16.17 


< 16.78 


17.70 ±0.01 


465 


152 


maybe 


790317 


17.59±0.19 


17.37±0.02 


16.94±0.11 


< 16.17 


< 16.78 


17.86±0.10 


573 


152 


maybe 


797384 


17.88 ±0.02 


17.47 ±0.03 


17.21 ±0.01 


16.23 ±0.08 


< 16.78 


18.09 ±0.01 


1109 


152 


yes 


799887 


17.60±0.18 


17.87±0.10 


17.50 ±0.06 


16.77 ±0.32 


17.27 ±0.42 


18.23 ±0.03 


817 


152 


maybe 


801865 


17.88 ±0.01 


17.66 ±0.06 


17.32 ±0.02 


16.38±0.10 


< 16.78 


18. 16 ±0.03 


698 


152 


yes 


803187 


18.17 ±0.03 


17.56 ±0.05 


17.28 ±0.04 


< 16.17 


< 16.78 


18.31 ±0.01 


1648 


152 


yes 


803471 


17.84 ±0.02 


17.58 ±0.03 


17.24 ±0.03 


16.41 ±0.04 


< 16.78 


18.11 ±0.01 


1253 


152 


yes 


806191 


17.92 ±0.04 


17.24 ±0.07 


17.12±0.05 


16.60±0.13 


< 16.78 


18.08 ±0.02 


657 


152 


maybe 


817663 


17.06 ±0.14 


17.28±0.13 


16.98±0.21 


< 16.17 


< 16.78 


17.54 ±0.08 


789 


152 


maybe 


Known Stars 


425399 


17.26 ±0.21 


< 16.79 


< 16.65 


< 16.17 


< 16.78 


17.38 ±3.38 


330 


152 


known stars 


564417 


17.68 ±0.63 


17.00 ±0.29 


17.01 ±0.40 


16.37±0.21 


< 16.78 


17.86±3.86 


1205 


152 


known stars 


619964 


17.37 ±0.32 


16.80 ±0.02 


16.80±0.17 


< 16.17 


< 16.78 


17.58 ±3.58 


1076 


152 


known stars 


660708 


17.18±0.40 


17.55 ±0.88 


17.23 ±0.68 


< 16.17 


17.41 ±0.63 


17.97 ±3.97 


1856 


151 


known stars 


696367 


17.72 ±0.66 


17.21 ±0.43 


16.91 ±0.27 


16.35±0.18 


< 16.78 


17.90 ±3.90 


499 


152 


known stars 



TABLE 5 

Gas-Phase Absorption Features of YSOs and possible YSOs 



SSTGC 




C 2 H 2 






HCN 






co 2 




C0 2 gas 


ID 




logWcol 


Abundance" 




logMol 


Abundance" 




logA'col 


Abundance" 


to solid ratio 




(K) 


(cm- 2 ) 




(K) 


(cm" 2 ) 




(K) 


(cm- 2 ) 






524665 


400 ±190 


16.9±0.1 


-5.4±0.1 


400 ±70 


17.0 ±0.2 


-5.3 ±0.2 


200 ±100 


17.2 ±0.3 


-5.1 ±0.3 


0.15 ±0.10 


726327 














100 ±50 


16.4±0.3 


-6.3 ±0.3 


0.04 ±0.03 


728480 


300 ±170 


15.7 ±0.3 


-6.9 ±0.3 








200 ±170 


16.2 ±0.2 


-6.4 ±0.2 


0.02 ±0.01 


761771 














100 ±50 


16.7 ±0.3 


-6.1 ±0.3 


0.02 ±0.01 


772151 














100 ±173 


16.6±0.3 


-6.2 ±0.3 


0.02 ±0.01 


797384 


100 ±160 


16.0 ±0.2 


-6.7 ±0.2 








100 ±50 


16.6±0.1 


-6.1 ±0.1 


0.03 ±0.01 


801865 








400 ±500 


16.5 ±0.4 


-6.3 ±0.4 


100 ± 158 


16.9 ±0.2 


-5.9 ±0.2 


0.05 ±0.03 


803187 


300 ±170 


16.4 ±0.2 


-6.4 ±0.2 


100±580 b 


16.4±0.6 b 


-6.4 ±0.6 


100 ±50 


16.8 ±0.2 


-6.0 ±0.2 


0.03 ±0.01 


803471 


200 ±160 


16.0 ±0.2 


-6.6 ±0.2 








100 ±150 


16.7 ±0.2 


-5.9 ±0.2 


0.04 ±0.02 



a Abundance relative to molecular hydrogren, log N/N(H 2 ). 

b Two local x 2 minima were found at T a = 100 K, \ogN„i = 16.4 and T a ss 700 K, logiV co i ss 16.7. 
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TABLE 6 

Ay FOR GC YSOS AND POSSIBLE YSOS 



SSTGC 


Ay (SL + LL) 


Ay (SH + LH) 


Ay (foreground) 11 


ID 


(mag) 


(mag) 


(mag) 


300758 


33.3 ±11.4 


33.7 ±23.6 


40±3 


304239 


28.7 ±7.0 


27.7 ±4.5 


43 ±4 


360559 


50.4 ±15.2 


48.5 ±9.9 


46±9 


370438 


53.3±6.8 


30.8 ±44.1 


62 ±17 


372630 


49.6 ±9.5 


36.0 ±4.1 


62 ±17 


496149 


38.8 ±7.4 


19.9 ±7.5 


31 ± 12 


524665 


20.5 ±10.7 


40.4 ±12.2 


43±9 


563780 


60.0 ±9.9 


39.7 ±3.3 


48 ±8 


610642 




24.6 ±9.5 


28±2 


618018 


34.3 ±9.6 


27.0 ±2.4 


27±5 


619522 


34.1 ±6.4 


48.5 ±3.3 


27±2 


653270 


17.6 ±0.4 


22.3 ±0.8 


22±3 


670953 


44.7 ± 1.5 


38.8±1.6 


31 ± 8 


679036 


51.3 ± 3.8 


43.3±4.1 


46 ±12 


718757 


31.9±8.4 


31.8±0.3 


76 ±22 


719445 


44.9 ±1.7 


53.5±2.1 


30±4 


722141 


39.5 ±1.9 


36.7 ±1.6 


19±5 


726327 


48.0 ±5.2 


43.3±2.1 


30±4 


728480 


44.3 ±4.0 


39.8 ±0.3 


39±11 


760679 


48.8 ±4.4 


48.9 ±0.8 


29±2 


761771 


69.5 ±7.3 


61.1 ±1.9 


47±8 


769305 


76.5 ±5.3 


53.5 ±2.3 


47 ±14 


770393 


30.4 ±1.0 


29.9 ±0.5 


36±5 


772151 


60.2 ±20.7 


64.5 ± 6.7 


24±3 


7 72981 


hoi I r» o 

48.3 ± 9.8 


40. 7 ± 0. 7 


36 ± 2 


773985 


70.6 ±6.9 


42.4 ±4.2 


36±2 


786009 


40.9 ±7.4 


34.2 ±0.4 


37 ±10 


790317 


45.4±14.9 


47.5 ±1.6 


31±4 


797384 


55.0±5.7 


55.1 ±0.9 


31 ± 1 


799887 


36.2 ±1.5 


36.2 ±2.3 


35 ±1 


801865 


66.7 ±9.4 


57.0 ±7.0 


31 ±4 


803187 


61.3 ±1.7 


57.6 ±1.4 


27±1 


803471 


46.8 ±2.1 


52.1 ±0.8 


31 ±4 


806191 


55.4 ±19.8 


55.7 ±0.9 


24±3 


817663 


32.6±14.0 


30.7 ±2.0 


29±1 


1 Based on the 2MASS and IRAC color-magnitude diagrams of GC red giant branch stars within 2' 


of the source (Schultheis et al.120091). 
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TABLE 7 

Photometric Properties of YSOs and possible YSOs 



SSTGC 




UKIDSS" 






IRAC b 




Synthetic 


SCUBA d 


ID 


J 


M 


K 


L3.6J 


[4.5J 


[5.8J 


L8.0J 


[24] c 


450 /^m 


850 fjxn 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(Jy) 


(Jy) 


O AA'V C O 

30075b 


i a f\f 1 A 1 1 

19.06 ±0.1 1 


1 A 1 A A 1 

14.77 ± 0.01 


1 O O I A A 1 

12.36 ± 0.01 


10.2 


9.0 


7.8 


6.2 


0.7 


12. 2± 1.4 


0.3±0.1 


304239 






16.88 ±0.18 


12.8 


10.2 


8.6 


7.7 


2.2 


< 5.0 


0.2±0.1 


360559 




IOTA 1 A OA 

Is. /4 ± 0.29 


1 A £Lf\ I A AT 

14.60 ± 0.03 


12.2 


11.1 


9.3 


7.6 


1 .7 






3 7043 o 






1 A f\0 1 A A A 

14.93 ± 0.04 


11.9 


10.0 


8.9 


7.4 


1 .8 






3 llbiU 




1 £. on 1 A AT 

16.39 ± 0.03 


1 T 1 A A 1 

13.63 ± 0.01 


10.3 


8.8 


7.7 


6.5 


1 .5 


1.9±0.8 


0.1 ±0.1 


ACiS" 1 A A 

496149 




1 AO 1 A 1 1 

16.93 ±0.11 


11 OA 1 A AO 

13.89 ±0.03 


12.0 


11.0 


9.6 


8.2 


1.3 






524665 






1 C 1 1 1 A 1 A 

15.71 ±0.10 


1 1 .4 


8.6 


7.1 


6. 1 


0.5 


428.5 ± 1.1 


36.9 ±0.1 


563 780 


1 O CjC 1 A 1 A 

18.56±0.10 


1 C Al 1 A A/1 

15.93 ± 0.04 


1 A AA I A A A 

14.09 ± 0.04 


1 1.6 


10.8 


8.2 




0.7 






610642 


1 O AT 1 A 1 C 

18.97±0.15 




1 O C ^ I A A 1 

12.56 ± 0.01 


9.8 


8.0 


6.6 


4.8 


—0.2 






618018 


1 C O O I A A 1 

15.33 ± 0.01 


1 /I CI 1 A A 1 

14.52 ± 0.01 


1 -1 r ~IC\ 1 A AO 

13. /9 ± 0.02 




9.6 


7.9 


6.5 


0.5 






619522 




1/1 C A 1 A A 1 

14.50 ± 0.01 


1 1 A"7 1 A A 1 

11.97 ± 0.01 


10.4 


9.3 


8.4 


7.7 


2.0 






6532/0 




1 £ <1 1 1 A AiC 

16.43 ± 0.06 


1 A AO I A A 1 

10.98 ± 0.01 


8.3 


7.1 


6.1 


5.7 


2.3 


26. 2± 1.6 


1.1 ±0.1 


670953 




1 T O O 1 A 1 C 

1 7.38 ± 0.15 


1 A AO I A f\£ 

14.92 ± 0.06 


1 1 .4 


8.9 


7.1 


6. 1 


0.6 






6 79036 


1 "7 AC I A m 

17.05 ± 0.02 


1 / f n I A A A 

16.50 ± 0.04 


1 A £.£. 1 A A A 

14.66 ± 0.04 


1 1.3 


9.4 


7.6 


6. 1 


0.7 


8.1 ±0.7 


1.3 ± 0.1 


/18 757 






1 C A 1 1 A A"7 

15.41 ± 0.07 


10.9 


9.3 


7.8 


6.0 


—0.3 






719445 




1 /" C A I A A A 

16.50 ± 0.04 


1 O C 1 1 A A 1 

13.51 ±0.01 


1 1.4 


9.1 


7.7 


6.0 


1 .3 






722141 






1 C O C I A A C 

15.35 ± 0.05 


13.0 


10.9 


9.3 


7.3 


0.3 






T~i£ "1 T7 

726327 


1 o o n I (i m 

18.29 ± 0.0/ 


11CC 1 A A 1 

13.55 ± 0.01 


11 CC I A A 1 

1 1.55 ± 0.01 


9.3 


7.9 


6.6 


4.9 








728480 






1 O OA I A A 1 

13.20 ± 0.01 


1 1.3 


10.3 


9.3 


7.5 


-0.3 


47.9 ±0.9 


8.0±0.1 


/606/9 






I r | A A "7 

15.6/ ±0.0/ 




10.3 


8.3 


6.5 


—0.2 


21.0±0.9 


0.9±0.1 


761771 






1 C A 1 1 A f\£ 

15.43 ± 0.06 


13.2 


10.6 




8.0 


1 .3 






"7^A7 Ac 

/69305 


18.50 ± 0.08 


1 O OO I A 1 A 

18.23 ±0.19 




11.8 


8.8 


7.2 


5.3 


—0.5 


492. 6± 1.2 


28.1 ±0.1 


/70393 




1 "7 1 "7 1 A ( \1 

1 1 A 1 ±0.0/ 


1 O O 1 1 A A 1 

13.81 ±0.01 


10.4 


8.7 


7.2 


5.1 








772151 


15.88 ± U.01 


14.85 ± 0.U1 


1 A n _1_ A AO 

14.22 ± 0.U2 


13.1 


11.7 


10.2 




2. 1 


568.3 ±1.4 


21.5 ±0.1 


772981 




16.11 ±0.03 


13.54±0.01 


11.2 


9.8 


8.2 


7.0 


0.0 






773985 






14.30 ±0.02 


12.0 


10.7 






-0.7 






786009 










10.7 






-0.6 


62.4 ±1.0 


2.0 ±0.1 


790317 




17.21 ±0.07 


14.40 ±0.02 




10.8 


9.2 


7.5 


0.5 


7214.4±3.5 


326.2 ±0.3 


797384 


18.52±0.08 


15.59 ±0.02 


13.72 ±0.01 




9.4 


7.7 


5.6 


-0.2 


340.4 ±3.6 


4.4 ±0.4 


799887 






14.24 ±0.02 


9.4 


7.2 


5.8 


5.3 


2.7 






801865 




17.07 ±0.07 






11.3 


10.3 


8.8 


3.3 






803187 


17.39±0.03 


16.61 ±0.05 


14.39 ±0.02 


12.2 


9.0 


7.2 


5.1 


-1.1 






803471 






13.50 ±0.01 


10.5 


8.8 


7.5 


6.0 


-0.1 






806191 


16.52±0.01 


15.43 ±0.01 


14.61 ±0.03 


12.6 


11.0 


9.1 


7.4 


1.6 


243.5 ±1.2 


1 1.3 ± 0.1 


817663 




16.35 ±0.03 


14.34 ±0.02 


12.6 


11.5 


10.1 


8.9 


2.3 


133.9±1.1 


8.1 ±0.1 



a Aperture3 magnitudes from UKIDSS DR2 i Warren et al. 20Q3). 

k Systematic errors of IRAC photometry were determined to be 0.1 mag, 0.1 mag, 0.15 mag, and 0.2 mags for channels 1, 2, 3, and 4, respectively, by comparing the Ra mirez et al.| 
1 2008) values with the measurements from the GLIMPSE II catalog i Churchwell et al. 2009). 
c Synthetic photometry based on IRS spectra. 

d Fluxes from SCUBA Legacy Catalogues I Di Francesco et al. 2008). 
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TABLE 8 

SED Fitting Results and Mass Estimates from Radio Continuum 



SSTGC 


logM* 


log Ltot 


A v (ISM) 


log M cm . 


log Age 


Radio Modeling 


ID 


V vl © ) 




(mag) 


(Mr^ vr~ l ) 


(vr) 


lr-irr AA t AA ^ \ 

lOgM* (Mq) 


References 


300758 


1 .03 ± 0.05 


3.64 ±0.23 


22.3 ±2.9 


-3.38 ±0.55 


4.43 ±0.47 


1.18 


i 


304239 


0.88 ±0.08 


3.18±0.28 


31.2 ±7.1 


-3.66 ±0.37 


4.44 ±0.81 






360559 


1.01 


3.37 


30.7 ±4.1 


-3.39 


4.12 






370438 


0.91 ±0.06 


3.22 ±0.27 


25.1 ±4.4 


-3.69 ±0.34 


3.89 ±0.78 






372630 


0.97 ±0.09 


3.34 ±0.24 


24.1 ±3.4 


-3.57 ±0.24 


3.29 ±0.22 






496149 


0.97 ±0.05 


3.51 ±0.19 


24.0 ±5.0 


-3.73 ±0.37 


4.72 ±0.75 






524665 


1.04 ±0.05 


3.90±0.16 


33.9±6.1 


-3.07 ±0.27 


4.94 ±0.31 






563780 


0.99 


3.45 


20.0 ±9.0 




4.53 






610642 


1.06 ±0.09 


3.84 ±0.29 


23.6±2.5 


-4.01 ±0.71 


4.95 ±0.96 






618018 


0.97 ±0.02 


3.35 ±0.10 


20.0 


-3.49 ±0.63 


4.53 ±0.20 






619522 


0.98 


3.40 


20.0 


-3.68 


4.58 






653270 


1.17 ±0.09 


4.32 ±0.24 


39.3 ±1.2 




6. 17 ±0.09 






670953 


1.37 ±0.14 


4.32±0.31 


33.8 ±7.7 


-2.72 ±0.27 


3.26 ±0.39 






679036 


0.98 ±0.01 


3.39 ±0.08 


20.0 


-3.74 ±0.55 


4.57 ±0.06 






718757 


1.13±0.07 


3.92±0.13 


33.1 ±5.8 


-3.52 ±0.45 


4.14±0.52 






719445 


0.88 ±0.07 


2.94 ±0.16 


22.0 ±2.3 


-4.19±0.18 


3.57 ±0.37 






722141 


1.04 ±0.05 


3.95±0.15 


35.8 ±5.0 


-3.22 ±0.27 


5.07 ±0.26 






726327 


1.36 ±0.05 


4.67 ±0.31 


20.6±1.1 


-2.85 ±0.45 


4.32 ±0.69 


1.27 


2 


728480 


1.18 ±0.05 


3.91 ±0.11 


22.7 ±3.8 


-3.47 ±0.20 


3.72 ±0.46 


1.26 


3 


760679 


1.42 


4.97 


40.0 




5.00 


1.26 


4 


761771 


1.00 


3.96 


20.0 




4.89 






769305 


1.24 ±0.01 


4.29 ±0.04 


20.0 


-3.04 ±0.09 


3.83±0.10 


1.27 


5 


770393 


1.25 ±0.06 


4.53±0.14 


35.9±2.8 


-4.08 ±0.11 


4.56 ±0.42 






772151 


1.14 


3.91 


20.0 




3.91 


1.26 


6 


772981 


1.20 ±0.03 


3.84±0.13 


20.0 ±0.1 


-3.39±0.18 


3.48±0.19 






773985 


1.21 ±0.09 


4.24 ±0.26 


21.8±1.9 


-3.17±0.27 


4.04 ±0.48 






786009 


1.15 ±0.06 


4.12±0.14 


33.9 ±6.9 


-3.58 ±0.36 


4.51 ±0.57 






790317 


1.13±0.07 


3.74 ±0.18 


21.7±2.0 


-3.40 ±0.25 


3.65 ±0.53 






797384 


1.20 ±0.04 


4.05 ±0.24 


20.0 


-3.07 ±0.34 


3.92 ±0.58 


1.29 


7 


799887 


1.11 ±0.04 


3.67 ±0.08 


20.0 


-2.24 ± 0.09 


3.19±0.06 






801865 


0.93 ±0.02 


3.15±0.18 


20.0 


-2.96 ±0.26 


4.70 ±0.09 






803187 


1.22 


4.26 


20.0 ±9.0 




3.80 


1.39 


8 


803471 


1.22 ±0.06 


3.93±0.16 


20.2 ±0.4 


-3.27 ±0.22 


3.49 ±0.48 






806191 


1.05 


3.61 


20.0 ±9.0 




4.17 






817663 


0.92 ±0.04 


3.01 ±0.21 


21.8±3.8 


-3.54 ±0.63 


4.34 ±0.36 







REF E RENCES. — Refe rences for radio observations: (1) Yusef-Zadeh et al. i 2009); (2) GPSR5 488-0 .028, Mehrineer et al. 1 1992); Becker end] 119941) ; (3) #8, Mehrin ger etaf] 
fi"992T); iMehringeJ fl99^>; (4) 2LC 000.563-0.044, [M ehringer et al. 1 1992); Yusef-Za deh et alj 12001; ILazio & Cordesl 120081); (5) GPSR 5 0.602-0.037, Becker et al l fl994l); 
Lazio & Cordes 1 1998); Yusef-Zadeh et al. 1 2004); Lazio & Cordes 1 2008); (6) 1LC 000.635-0.020, Lazio & Cordes 1 1998); White et al. i 2005); (7) SGR B2 HII P, Mehringer et al. 
1 1993); Yusef-Zadeh et al. 1 2004); (8) GPSR5 0.693-0.046, Zoonematkermani et al. (1990); Becker et al. 1 1994); Lazio & Cordes 1 1998); Yusef-Zadeh et al. 1 2004); White et al. ( 2005); 
Lazio & Cordes 1 2008). 



